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PREFACE 



In recent years great progress has been made in the gas turbine 
field. The improvements made in high temperature materials, compressor 
and turbine design, and compact heat exchangers have made the gas tur- 
bine power plant a serious competitor of steam plants and diesel engines 
for many applications. There are, however, many problems still exist- 
ing that tend to make conventional gas turbines unattractive. 

The use of residual fuels in the open cycle results in ash deposits 
and corrosion in the turbines which seriously limits the operating life 
of the units. Furthermore, open cycle plants are limited in power rat- 
ing since the components become too large for outputs greater than about 
20,000 HP. The closed cycle gas turbine plant eliminates these problems 
since only clean air circulates through the cycle and high pressures 
permit components of a practical size for outputs much larger than 
20,000 HP. For a ship propulsion unit where the power requirements 
vary over a wide range, the closed cycle has the additional advantage 
of good part load efficiency. A major factor in the slow progress of 
the closed cycle plant has been the inability to manufacture a reliable 
air heater cheaply enough to permit the gas turbine plant to compete 
with conventional steam plants. 

The future of the gas turbines may depend on how well they can be 
adapted for use with nuclear power. At the present time the only ship 
propulsion plant employing a nuclear reactor is essentially a steam 
plant. The size of this plant is large since a two loop system is 
required in order that energy may be conveyed from the reactor to the 



ii 



secondary working fluid. A simpler and smaller plant would be one in 
which there is only one fluid serving both as the working fluid and 
reactor coolant. Such a plant is a closed cycle gas turbine plant 
using helium at high pressure. Helium has good coolant characteristics 
including a negligible cross section for neutron absorption. If pure, 
it is not subject to induced activity. By proper design of the cooling 
channels and filtering system, fission products can be prevented from 
being carried through that portion of the cycle which is external to 
the reactor. This would require only the reactor to be shielded. It 
has been presumed up to this point that helium is a satisfactory gas 
turbine fluid. Although helium turbines have never been built, the 
following general statements can be made by comparing the properties of 
helium with air: 

(1) Better heat transfer characteristics of helium results in 
smaller heat exchangers. 

(2) Lower molecular weight of helium requires more stages in the 
turbo-machines for the same peripheral speed. 

(3) Higher velocity of sound in helium permits higher peripheral 
speeds before Mach effect occurs. 

It would appear from these statements that the advantage of smaller 
heat exchangers for a helium cycle may be more than offset by requiring 
more expensive turbines and compressors as well as a costlier and less 
available gas. Design criteria for helium compressors are based on 
stress limitations, whereas for air compressors the Mach effect becomes 
the limiting factor. Thus, the higher peripheral speed permitted for 
the helium compressors will partially counteract the effect of lower 
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molecular weight on the number of stages required. Furthermore, improv- 
ed efficiency and fewer stages of compression may result by employing 
a higher effectiveness regenerator than would be practical with air. 

It is the writer's objective to investigate a 50,000 HP closed- 
cycle helium gas turbine power plant for ship propulsion and determine 
the approximate sizes of the major components. Unfortunately, the 
scope of this report does not cover many of the interesting problems 
which would have to be investigated and solved before this plant could 
become a reality. Some of these are discussed briefly in Chapter IX. 

It has been necessary to include a list of symbols at the begin- 
ning of each chapter since the same symbols may have different meanings 
in different chapters. 

The writer wishes to express his deep appreciation for the encour- 
agement and assistance given him by Professor M. H. Vavra and Professor 
H. M. Wright in the preparation of this report. The writer is also 
grateful for the fine clerical support given him by his typist, Mrs. 
Marilyn Heard. 



iv 




' 



TABLE OF CONTENTS 



Item Title Page 

Chapter I Thermodynamic Properties of Helium .... 1 

Chapter II Preliminary Cycle Analysis 12 

Chapter III Heat Exchangers 25 

Chapter IV Compressors 51 

Chapter V Turbines 82 

Chapter VI Reactor 105 

Chapter VII Installation 115 

Chapter VIII Stresses 133 

Chapter IX Summary and Conclusions 160 

Bibliography 163 



v 



LIST OF ILLUSTRATIONS 



Figure Title Page 

1-1 Variation of Absolute Viscosity and Thermal 

Conductivity with Temperature 7 

1-2 Variation of B q with Temperature 8 

1-3 Variation of Enthalpy with Pressure for Given 

Temperatures 9, 10 

1- 4 Variation of ~\j/ ^ and y/ with Initial 

Pressure 11 

2- 1 Schematic Diagram of a Closed Cycle Gas 

Turbine Plant 15 

2-2 Effect of Pressure Ratio and Regenerator 

Effectiveness on Plant Thermal Efficiency . . . 22 

2-3 Effect of Pressure Ratio on Helium Rate .... 23 

2- 4 Effect of Pressure Ratio on Work Ratio .... 24 

3- 1 Sketch of a Bar Fin 31 

3-2 Sketch Showing Turbulent Boundary Layer 

between Tube Fins 36 

3-3 Heat Exchanger Tube 37 

3-4 Tube Arrangement 37 

3-5 Schematic Representation of Regenerator .... 40 

3-6 Schematic Representation of Precooler 42 

3-7 Schematic Representation of Intercooler .... 43 

3-8 Variation of Regenerator Characteristics with 

Average Helium Velocity 47 

3-9 Variation of Precooler Characteristics with 

Average Helium Velocity 48 

3-10 Variation of Intercooler Characteristics with 

Average Helium Velocity 49 

vi 


















' 



< 



Figure Title Page 

4-1 Schematic Representation of an Axial Flow 

Compressor 54 

4-2 Typical h-s Diagram for an Axial Flow 

Compressor 55 

4-3 Schematic Representation of an Axial Flow 

Compressor Stage and a Corresponding Velocity 
Diagram 57 

4-4 Typical h-s Diagram for an Axial Flow 

Compressor Stage 58 

4-5 Dimensionless Velocity Diagram for an Axial 

Flow Compressor Stage 59 

4-6 Nominal Values of Drag Profile Coefficient . . 62 

4-7 Nominal Values of Fluid Deflection 62 

4-8 Nominal Values of Fluid Deflection 62 



4-9 T-S Diagram Illustrating Reheat Factor • . . . 67 



4-10 

4-11 

4-12 



Effect of (3 2 on Stage Efficiency for Given 

Values of (s/c) 70 

R 

Effect of @2 on Rotational Speed for Given 

Values of (s/c)„ and h 72 

' R n 

Effect of U on Number of Stages and Last 

Stage Blade Heights 78 



5-1 



Schematic Representation of an Axial Flow 

Turbine 84 



5-2 Typical h-s Diagram for an Axial Flow Turbine 86 

5-3 Schematic Representation of an Axial Flow 

Turbine Stage and Corresponding Velocity 
Diagram 88 

5-4 Typical h-s Diagram for an Axial Flow Turbine 

Stage 88 

5-5 Typical h-s Diagram for Relative Flow in 

Axial Flow Turbine Rotor Blades 90 



5-6 



Velocity Diagrams for 507. Reaction and Modified * 
Rateau Stages 95 



vii 



Figure Title Page 

5-7 Effect of on Turbine Stage Efficiency . . 97 

2 

5-8 Effect of on Number of Turbine Stages 

Required 2 98 

7-1 Channel of Constant Flow Area 117 

7-2 Channels of Decreasing Flow Area 117 

7-3 Channel with an Abrupt Increase in Flow Area . 118 

7-4 Diffuser 118 

7-5 Regenerator and Precooler 122 

7-6 Schematic Representation of Flow through the 

Regenerator and Precooler 123 

7-7 Schematic Representation of Flow through the 

Regenerator 123 

7-8 Intercooler 126 

7-9 Schematic Representation of Flow through the 

Intercooler 126 

7-10 Reactor 129 

7-11 Schematic Representation of Flow through the 

Reactor 129 

7- 12 Power Plant Arrangement 132 

8- 1 Forces Acting on a Compressor Rotor Blade . . . 135 

8-2 First Stage Rotor Blade for Compressor No. 1 . 140 

8-3 First Stage Stator Blade for Compressor No. 1 . 141 

8-4 Momentum Theorem Control Surfaces for the 

Blades of Compressor No. 1 143 

8-5 Forces Acting on a Turbine Rotor Blade .... 147 

8-6 First Stage Rotor Blade Profiles for 

Turbine No. 1 150 

8-7 Momentum Theorem Control Surface for a 

Turbine Rotor Blade 153 



viii 



Figure Title Page 

8-8 Turbine Nozzle Profile 153 

8-9 Turbine Disk 155 

8-10 Approximate Cross Section of Turbine Disk Rim . 156 

8-11 First and Second Stages of Turbine No. 1 . . . . 159 



ix 



LIST OF TABLES 



Table 

2-1 

2-2 

3-1 

3-2 

3-3 

3-4 

3-5 

3-6 

3-7 

3-8 

3-9 

3- 10 

4- 1 
4-2 

4-3 

4-4 

4-5 



Title 



Plant Performance as a Function of Pressure 



Ratio 20 

Thermodynamic Properties for Cycle Having a 
Pressure Ratio of 2.25 and a Regenerator 
Effectiveness of 0.9 21 

Properties of Helium Entering and Leaving 
Regenerator 40 

Average Properties of Fluids in the 

Regenerator 41 

Properties of Helium Entering and Leaving 
Precooler 42 



Average Properties of Fluids in the Precooler . 42 

Property of Fluids Entering and Leaving 



Intercooler 43 

Average Properties of Fluids in the 

Intercooler 43 

Tube Dimensions 45 



Heat Exchanger Areas and Associated Quantities . 45 

List of Constants for the Regenerator, 

Precooler, and Intercooler 46 

Data for Selected Heat Exchangers 50 

Effect of (s/c)^ on Compressor Blading 71 

Data for First and Last Stages of Compressor 

No. 1 73 

First Approximation Data for First and Last 
Stages of Compressor No. 1 74 

Additional Data for Compressor No. 1 75 

Final Data for First and Last Stages of 
Compressor No. 1 77 



x 



Table Title Page 

4-6 First Approximation Data for First and Last 

Stages of Compressor No, 2 79 

4-7 Additional Data for Compressor No. 2 80 

4- 8 Final Data for First and Last Stages of 

Compressor No. 2 81 

5- 1 First and Last Stage Data for Turbine No. 1 . . 101 

5- 2 First and Last Stage Data for Turbine No. 2 . . 102 

6- 1 Summary of Reactor Core Data 114 

7- 1 Piping Data 120 

7-2 Summary of Pressure Losses in the Regenerator 

and Precooler 124 

7-3 Summary of Pressure Losses in the Regenerator 

(H. P. Side) 125 

7-4 Suiranary of Pressure Losses in the Intercooler . 127 

7-5 Summary of Pressure Drops in the Reactor .... 128 

7-6 Comparison of Assumed and Calculated Values 

of AJi 130 

P. 

i 

7- 7 Comparison of Assumed and Calculated Values 

of Machine Efficiencies 131 

8- 1 Velocity Diagram Data for the First Stage of 

Compressor No. 1 137 

8-2 Velocity Diagram Data for the Last Stage of 

Compressor No. 1 138 

8-3 Airfoil Section Properties for a Chord of 

0.82 Inches 139 

8-4 Values of 'T and 5 139 

8-5 Centrifugal Forces and Corresponding Stresses 

for the Rotor Blades of Compressor No. 1 ... . 142 

8-6 Gas Forces Acting on the Blades of 

Compressor No. 1 145 

xi 



Table Title Page 

8-7 Blade Root Stresses Due to Gas Forces 146 

8-8 Velocity Diagram Data for First Stage of 

Turbine No. 1 148 

8-9 Centrifugal Loading on the Turbine Disks . . . 157 

8- 10 Widths of Constant Stress Disk 158 

9- 1 Summary of Descriptive Data for the Major 

Components 161 



xii 



CHAPTER I 



THERMODYNAMIC PROPERTIES OF HELIUM 



Symbols used in this chapter: 

Symbol Description 

B virial coefficient 

o 

Cp specific heat at constant pressure 

c v specific heat at constant volume 

f reheat factor for oo stages 

h enthalpy 

k thermal conductivity 

P pressure 

R gas constant 



T 

t 

v 

r 

A 

V 

H 

r 



absolute temperature 
temperature 
specific volume 
ratio of 

c 

v 

increment 

efficiency 

absolute viscosity 

corresponds to Y -1 
T 



Subscript Definition 
o reference quantity 

c compressor or compression process 



Units 
ft 3 /lb 
BTU/lb - °F 
BTU/lb - *P 



BTU/lb 

BTU/ft 2 - hr - ®F/ft 

psia 

ft - lb f 
lb 

m 

*R 

e F 

ft 3 /lb 



lb/hr - ft 



1 



Subscript Definition 



i initial 

f final 

s isentropic process 

t turbine or expansion 

General 

In studying a cycle that uses fluids such as steam or air, or one 
that behaves like an ideal gas, there are tables, charts and equations 
which can be used to obtain accurate thermodynamic information. This 
is not the case for helium. Several equations of state have been 
proposed based on theory and empirical data obtained at low tempera- 
tures. Which of these equations, if any, is applicable at higher 
temperatures will be known only after more extensive empirical data 
are available. Until then, extrapolation from the lower temperatures 
will have to suffice. The extrapolation methods used in this investi- 
gation are based on the research done by Akin ^ \ His tables for 
enthalpy, entropy, and specific volume extend to 6Q0°F whereas the 
highest temperature in the cycle being investigated is 1350°F. His 
extrapolation formulas for absolute viscosity and thermal conductivity 
extend to temperatures beyond 135Q°F. 

Absolute Viscosity , yU 

(2 ) 

Akin found Keesom's ' 7 equation for yU accurate within 1% for 
temperatures between -450 and 2000°F. This equation, plotted in 
Figure 1-1, is: 

0.647 /t 

^ = d.3 1 5 T x 10 tb/hr-ff ( 1 - 1 ) 
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Thermal Conduct ivity , k 

The following equation is based on kinetic theory and test results 
at temperatures below 212°F. It is plotted in Figure 1-1. 

K = A< Cy € QTU/if-hr - °F/ff 

or K ~ I-8I6& /u -for t > 80 °F (1-2) 

Specific Volume , v 

Akin 's, results on specific volume are based on the Beattie-Bridge- 
man equation of state. Comparisons with the co-volume equation of 
state show good agreement in the range of pressures and temperatures 
encountered in this report. Because of its simpler form, the latter 
equation has been used throughout. 

P = P S,Q 0-3) 

- T - 

where ft, = 0. 306 ) - 0.822 (l-3a) 

3 

and P = 2,6 829 xtr lj>£— 

I b») ~ °p - ir\ 2 

Equation (l-3a) is plotted in Figure 1-2. The specific volume for 

temperatures from 0°F to 1400°F when calculated from equation (1-3) 

3 

is 0.047 to 0.040 ft /lb greater than that obtained by assuming ideal 
gas behavior at the same temperature. 

Enthalpy , h 

For enthalpies below 600°F, Akin’s tables have been used. The 
extrapolation formula for a constant pressure process can be expressed 
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as : 



h t = k + c p 0 (t 'to) + f(± j P) 

where h 0 = h ot ~ 44 0 F 
and Cp' = /. 2 4/ B T {jy / 1 b - F 



(1-4) 



At temperatures above 400°F it is apparent from the tables that for a 
constant pressure process: 



Figure 1-3 contains plots of h vs, P for temperatures from 100°F to 
600°F, 

Isentropic Work of Expansion and Compression 

For an ideal gas the work of compression for a reversible-adiaba- 
tic process is: 




where > 400 F 

and C p =U+ BTU/th-'F 



Since extrapolation is necessary above 600°F # (1-3) becomes: 



K = + /.2+(t-600) BTCj/lb a-5.) 



~f~ I 




( 1 - 6 ) 



4 



Similarly, the work of expansion is: 



TZ c p 




BTU/ lb 



(1-7) 



An analysis of Akin's tables reveals that these equations are not 



applicable for helium if both c and y are assumed to be constant. 

P 

Since c^ is constant (1.24 BTU/lb °F) equations (1-6) and (1-7) were 

•y -1 p 

solved for — — using various combinations of T^, P^, and f/P^. 



V-l , 



The results indicated that ~~ is a function of P, only. Since this 



function differs for the compression and expansion processes, 



T -1 



has been replaced by y/ c and y' respectively. 



* h s = 



y 

I) - 



BTU/lb 



<l-6a) 



A f) s = T l 



I - % 



ff\ Yt 



{_ 

Pi 



BTU/lb 



(l-7a) 



Figure 1-4 contains plots of and yy vs. P^. 

t 

Reheat Factor . 1+foo 

The values of ys and y/ ^ are also applicable to expressions for 

evaluating reheat factors. Assuming ideal gas behavior and replacing 
Y -1 

by y or y/ , the reheat factor, 1+foo # becomes: 

Y c t 

for compression 

( 1 - 8 ) 




5 



1+f 



oo 



for expansion 



% 



I -[%. 



g}M 

Pa 



I - 



w 



(1-9) 



A more detailed analysis of reheat factors is given in Chapter IV on 
"Compressors” . 

By means of the aforementioned formulas and/or the plots, all the 
thermodynamic properties that are required in the following cycle 
analysis can be obtained. 
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Thermal Conductivity, k, (BTU/hr - ft z - °F/ft) x 10 




<20X20 PER INCH ’*\0R 





CHAPTER II 



PRELIMINARY CYCLE ANALYSIS 



Symbols used in this chapter: 
Symbol Description 

B virial coefficient 



o 



c 

P 



c 

V 



A H 



h 

k 

• 

m 

P 

R 

T 

t 

v 

V 



specific heat at, constant pressure 

specific heat at constant volume 

total isentropic enthalpy change for 
the compressor or turbine 

enthalpy 

thermal conductivity 
mass rate of flow 
pressure 
gas constant 

absolute temperature 
temperature 
specific volume 
ratio of 

c 



v 



Units 
ft 3 /ib 
BTU/lb - °F 
BTU/lb - °F 

BTU/lb 

BTU/lb 

BTU/ft 2 - hr - 

lb/sec 

psia 

ft - lb f 

lb ^~°F 
m 

°R 

# F 

ft 3 /lb 



4 

V 

r 



increment 
efficiency 
absolute viscosity 
corresponds to — ; 



lb/hr - ft 



°F/f t 
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Subscript Definition 



o 



reference quantity 



c 



compressor or compression process 



i 



initial 



f 



final 



R 



regenerator 



s 



isentropic process 



t 



turbine or expansion 



Th 



thermal 



General 



Since it is the objective of this report to analyze a 50,000 HP 
closed-cycle gas turbine plant, it is of prime importance to select a 
set of conditions which will result in a plant of high thermal effi- 
ciency and reasonable size which is still practical and economical to 
build, operate and maintain. To make such a selection one must know 
the following: 

(a) Components of the cycle 

(b) Reasonable estimate of the machine efficiencies 

(c) Anticipated losses associated with the various components and 
interconnecting piping. 

(d) Regenerator effectiveness 

(e) Limiting temperatures, both upper and lower 

(f) Highest pressure permissible in the cycle. 

Components of the Cycle 

This is not a study of various types and arrangements of cycles, 
hence a cycle of minimum complexity was chosen. The basic cycle is 
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( 3 ) 



similar to the conventional plant built and operated by Escher-Wyss 
and consists of a single heat source, regenerator, precooler, compress- 
ors with intercooling, and two turbines. For the heat source, the 
conventional air heater is replaced by a homeogeneous reactor, the core 
consisting of a graphite and highly enriched uranium mixture. Since 
auxiliary equipment, including control apparatus, has little effect on 
the thermodynamics of the cycle it is not included here. A diagram- 
matic sketch of the basic cycle is shown in Figure 2-1. 

Machine Efficiencies 

Assigning efficiencies to the two compressors and two turbines 
can at best be only a first approximation. Turbo -machines are built 
today with overall efficiencies over 80% when used in conventional 
cycles. However, for a given cycle such as the one being investigated, 
practical design considerations may preclude these high efficiencies 
and a conservative estimate is in order. No distinction is made in 
the preliminary analysis between internal and overall machine effi- 
ciencies and no separate mechanical efficiency is specified. The 
efficiencies for Compressor No. 2 and Turbine No. 2 are lower than 
for the No. 1 units in the event that radial f low turbo-machines are 
practical . 



MACHINE 



EFFICIENCY (Assumed) 



Turbine No. 1 



0.85 



Turbine No. 2 



0.80 



Compressor No. 1 



0.84 



Compressor No. 2 
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0.80 



c*o 
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Figure 2-1 Schematic Diagram of a Closed Cycle Gas Turbine Plant 





Pressure Drops 



The losses associated with each component are expressed here as a 
rsrtio of the static pressure loss across the component to the inlet static 
pressure. Again, these are only estimates based on rough calculations 
and experience and include not only the drop across the component but 
also associated inlet and exit piping effects. 

COMPONENT (Assumed) 

i 



Reactor 0.010 

Regenerator 



High Pressure Side 


0.015 


Low Pressure Side 


0.010 


Precooler 


0.0075 


Intercooler 


0.005 


Regenerator Effectiveness 





The regenerator effectiveness is defined as: 




t±5 zh± 

h 8 -h+ 



( 2 - 1 ) 



A high effectiveness is usually desirable from the standpoint of high 
thermal efficiency. It is not practically feasible to have tj greater 
than 0.9 because of the large amount of surface required. la the 
preliminary analysis is taken as 0.7, 0.8 and 0.9 to indicate its 
effect on efficiency. 

Limiting Temperatures of Working Fluid 

The lower temperature is dependent primarily on the cooling medium, 
in this case sea water. For a given plant a decrease in the lower 
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temperature will improve plant efficiency. However, since a ship may 
operate in tropic waters as well as arctic waters, a conservative 
lower limit of 100°F is chosen. The upper temperature limit depends 
on materials and is gradually increasing as materials and cooling 
techniques are improved. Temperatures of 1450 and 1500°F are common 
in present-day gas turbine cycles, but because of the high pressure 
encountered in this closed cycle an upper temperature of 1350°F is 
assumed , 

Maximum Pressure 

The maximum pressure in the cycle has a slight effect on effi- 

> 

ciency due to the variation of anc * pressure. More 

important, however, is the effect of pressure on the size of the plant. 
Higher pressures mean smaller specific volumes, hence smaller flow 
areas. Based on steam plant design a maximum cycle pressure of 1000 
psia seems reasonable for a plant of this size. 

Cycle Analysis 

Now, with the foregoing assumptions, together with the thermody- 
namic information from Chapter I, the cycle of Figure 2-1 has been 
analyzed by standard methods to obtain thermodynamic properties at 
each point for pressure ratios between 1.25 and 5.0. In addition, 
the following performance indices are plotted in Figures 2-2, 2-3 and 
2-4 respectively: 

(a) Thermal efficiency based on work output at turbine shaft 

(b) Helium rate 

(c) Work ratio 
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A list of equations used in the analysis follows: 



AH = 
aH = 
h t = f 
aH s = 
AH — 

%= 

Helium 



• 68 ^- r + 6o f t 3 /ib 

AHi 

% 

AH s * t) t 



>6o o-F + A2f (t -600) BTU/lb 



X 



/.2 4 T; 



BTU/lb 



> f 



(.24 T 



t 



BTU/lb 






Rate - 



2545 

^7’^t 



/jb/ hp-h> 



(1-3) 



( 2 - 2 ) 



(2-3) 



<l-5a) 



(l-6a ) 



(l-7a) 



(2-4) 



(2-5) 
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Work Ratio 



( 2 - 6 ) 



hi zJl 

h G -h 



8 



n — b s-iJiA 

»R ~ L, L 

^8 h 4 - 



(2-1) 



Conclusion 

In studying Figure 2-2 we observe the following points of maximum 
efficiency: 



*2_r 


P 

4/P l 


' ?.h 


0.7 


3.00 


0.268 


0.8 


2.65 


0.292 


0.9 


2.25 


0.329 



If this cycle is to compete with other power plants it must have a 
higher thermal efficiency and/or be smaller in size. A high effective- 
ness regenerator provides the higher efficiency at the expense of 
increased size of plant. At the lower effectiveness, not only is the 
thermal efficiency less, but the larger pressure ratio will result in 
several more stages for the turbo-machines. Placing emphasis on the 
higher efficiency and fewer stages, the remainder of this report consid- 
ers only the cycle having » 2.25 and r 3 0.9. Table 2-2 

lists the thermodynamic properties at each point in this cycle. 
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TABLE 2-1 



PLANT PERFORMANCE AS A FUNCTION OF PRESSURE RATIO 



Pressure 

Ratio 

P 

r 4/P t 


Plant Thermal Efficiency, y th 
= 0.9 R = 0.8 R = 0.7 


Helium 

Rate 

( Ib/hp-hr ) 


Work 

Ratio 


1.25 


0.184 


0.122 


0.091 


51.7 


0.378 


1.50 


0.279 


0.210 


0.169 


24.4 


0.413 


1.75 


0.312 


* 0.253 


0.212 


17.7 


0.410 


2.00 


0.325 


0.274 


0.237 


14.8 


0.400 


2.25 


0.329 


0.285 


0.252 


13.2 


0.388 


2.50 


0.327 


0.290 


0.261 


12.2 


0.375 


2.75 


0.323 


0.291 


0.265 


11.6 


0.362 


3.00 


0.318 


0.291 


0.268 


11.1 


0.351 


4.00 


0.290 


0.275 


0.261 


10.4 


0.308 


5.00 


0.261 


0.252 


0.244 


10.6 


0.270 



2.0 



TABLE 2-2 



THERMODYNAMIC PROPERTIES FOR CYCLE HAVING *4^ = 2.25 AND Kj R = 0.9 



Location 
(Fig. 2-1) 


P (psia) 


h(BTU/lb) 


t(°F) 


v(ft /lb) 


1 


444 .4 


701.5 


100 


3.428 


2 


666.7 


848.8 


216.9 


2.768 


3 


663.4 


703.7 


100 


2.309 


4 


1000 


861.9 


224.9 


1.883 


5 


985 


1667.8 


874.8 


3.677 


6 


975.2 


2256.9 


1350 


5.019 


7 


631.2 


1951.0 


1106.3 


6.692 


8 


452.2 


1757.3 


951.3 


8.413 


9 


447.7 


951.4 


301.5 


4.607 


m = 


182.4 


lb / sec 


(For Net output 


= 50,000 HP) 


^ th = 


0.329 








Helium 
Rate * 


13.2 


lb /hp-hr 






Work 

Ratio 


0.388 
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CHAPTER III 



HEAT EXCHANGERS 



Symbols used in this chapter: 



Symbol 


Description 


Units 




A 


heat transfer area per tube 


ft 2 




a 


free flow area per tube 


ft 2 




a t 


frontal area per tube 


CM 

U 

U-l 




a 


fin longitudinal length 


in. 




c 


fin thickness 


in. 




c 

p 


specific heat at constant pressure 


BTU/lb - 


°F 


D 


heat exchanger shell inside diameter 


ft 




D h 


hydraulic diameter 


ft 




d 


tube diameter 


in. 




h 


convection heat transfer coefficient 


BTU/ft 2 


- hr 


k 


thermal conductivity 


BTU/ft 2 


- hr 


L 


heat exchanger tube length 


ft 




X 


fin radial length 


in. 




m 


(defined by equation 3-10a, page 33) 


— 




m 


mass rate of flow 


lb/sec 




N 


total number of tubes 


— 




Nu 


Nusselt number 


— 




P 


pressure 


psia 




Pr 


Prandtl number 


— 




Re 


Reynolds number 


— 




s 


tube center to center spacing 


in. 





°F 

°F/f t 



25 



Symbol 


Description 


Units 


t 


temperature 


°F 


U 


overall heat transfer coefficient 


BTU/f t 


V 


velocity 


ft /sec 


V 


specific volume 


ft 3 /lb 


w 

6 


. . (me) 

ratio : p c ' 

(m c ), 
p h 


— 


v;all thickness or boundary length 
thickness 


in . 


A 


increment 


— 


t 


friction factor 


— 


% 


fin efficiency 


— 


% 


fin temperature effectiveness 


— 


p 


absolute viscosity 


Ib/hr 


p 


density 


lb/ft 3 


<p 


dimensionless overall heat transfer 
coefficient 


— 


y 


dimensionless convection heat 
transfer coefficient 


— 


SI 


A 

ratio of heat transfer areas, c 


— 


5 


ratio of fluid flow areas, 


— 



Subscript 


Definition 


c 


cold fluid 


e 


exit 


f 


fin 


h 


hot fluid 
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Subscript Definition 



i 



inlet, inside, or initial 



o 



outside 



t 



tube 



w 



wall 



General 



Having selected a pressure ratio of 2*25 from the preliminary 
analysis, the components will now be more closely investigated to 
determine if they can be of practical design. Since one of the more 
important requirements is the size of the plant, a design study of the 
regenerator, intercooler and precooier is necessary. Before going into 
this study, some review of heat transfer theory as used herein may be 
advantageous • 

Heat Transfer Theory 

In the usual type of heat exchanger heat transfer occurs between 
two* fluids of different temperatures separated by a wall of finite 
thickness. The quantity of heat, dQ, transmitted through an element of 
area, dA, can be expressed as: 



where U is the overall heat transfer coefficient associated with are* 
dA. 

If the wall surface is a plane, i.e., no extended surface: 




(3-1) 




(3-2) 
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A » ^ 



If U = U 



h 1 



U L 



i_ + s. 






Ac K 



(3-2a) 



Ah 



.( 4 ) 



Kays v y has modified (3-2a) to include extended surfaces on both 
sides of the wall: 



and 



I _ I 



Ul 



^°>A 



Ay) 

A, 



a; m 



where = effective surface area on hot side 



Oo A * effective surface area on cold side, 
' c c 



(3-3) 



For thin walled surfaces the terra 



w 






is negligible: 



_ I 



~A^°^ C 



(3-3a) 



U h *l°h\ 

Combining the techniques of Kays and Traupel v "' / it is convenient 
to change (3-3a) to a dimensionless form: 



(5) 



<P h \\ WcW * 1 



(3-4) 



where 



t 



u 

l° c pV 



(3 -4a) 



28 



* 



y 



h 

(°^V 



W = 



(^ C p)c 

c f\ 



(also known as 
Stanton Number) (3-4b) 



(3-4c) 



5 = ( J-4d) 

U h 

O. — (3-4e) 

A h 



Considering the hot fluid only, the heat loss, dQ, for a tempera- 
ture change dt can be expressed as: 

dQ h = = -(pQVcJc U h (3-5) 



Combining (3-5) with (3-1) 

dAh __ 0 Ith 



<t >, I 



h a 



L — _ 
h 






(3-6) 



If we assume that the configuration and dimensions of the heat 

exchanger elements do not change along the direction of the flow, then 

dA^ e dx, where x is a lengthwise coordinate, L is length of the 

L 

surface and is the total hot side area. Integrating equation (3-6): 




dt 



(3-7) 



o 



or 



A, 

O h 





_i 

V*c 




(3-7a) 



(p^ will change as the fluid flow proceeds along the surface and 
its function of x must be known before the integral can be evaluated. 
Likewise, t^ - t^ changes and its function of t^ must be known. 
Usually these functions are not known and it is common practice to 
simplify (3-7) as follows: 



Al A = (U-te\ 
O h ^ At 



(3-7b) 



where 




^d(i) 



so 

and == — 

/ A \ 

Bowman, Mueller and Nagle v ) have determined relations for At 
for many types of heat exchangers. For the longitudinal counter-flow 
exchanger A t is the familiar log mean temperature difference: 



ti-u 



I 



At 

Wjt C 



^ _ (bu -tc e 1 - 

In (t ht -Q 

(k-Q 



(3-8) 



Close approximation is obtained if (p is evaluated by equation 
(3-4) using average values of ~\f/ and -y/ ^ . 
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Up to this point the type of extended surface baa not been speci- 
fied. Once this is specified there remains the problem of evaluating 



>lo ■ 

h A has been defined as the effective heat transfer area. For a 
/° 

finned surface this effective area can be expressed as the sum of the 
wall area (not including fin base), and an effective fin area. 

From the definition of fin efficiency, IJ this effective fin area is 
7f ^f w ^ ere Af is the total surface area of the fin. 

fyA = Aw 4- ry A, 

Aw — A — Af 

•’* % ~ 1 ~ ^ (3-9) 



For several types of fins the solutions for rj ^ are known. 
(Gardner^ 

As an example, for a bar fin of constant thickness shown in 
Figure 3-1: 




7 f 



ta_n_h (jjnj) 
mX 



(3-10) 



where tf/~\ — U ?JxCLJi£l 
> Af L C 
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if L » c : 



rn 



(3-10a ) 




It is noted that Gardner’s solutions for rj ^ are based on several 
assumptions. Of these, the following may result in appreciable error: 

(1) The heat transfer coefficient, h, is uniform over the 
entire fin surface. 

(2) The temperature of the base of the fin is uniform. 

(3) Temperature gradients normal to the fin surface are 
neglected. 

(4) Heat transfer through the outermost end is negligible. 

The only type of heat exchanger involved in this investigation is 

tubular with external fins. By introducing the hydraulic diameter, D^: 




40i 

K/l 



equation (3-7b) becomes: 



L 


1 


(ti t e \ 


D* ' 


’ +% 


At 



(3-11) 



(3-12) 



The above equation will be one of the principal equations employed 
in this investigation for design of the heat exchangers. 

Pressure Drops 

Fluid pressure drop in the heat exchanger must be considered in 
addition to heat transfer aspects. It has been found both experimen- 
tally and analytically that the two are interdependent. Improvement in 



32 



heat transfer characteristics results in increased pressure drop, so a 
compromise has to be made in most designs. The relation used by 
Traupel^\ McAdams^^ and others for longitudinal flow inside of tubes 
is : 



AP=! -L _vl 

A 2 9 v 



(3-13) 



where £ is the coefficient of resistance. It is identical to ^-f 
in McAdams . 

For turbulent flow through smooth tubes'" 



r 



_ ” O.Z 

0.18 4 Re (3-14) 

Sooo < Re < SOQjOoo 



Similarly, for determining pressure drop outside of the tubes (on 
finned side), equation (3-13) and (3-14) are applicable if the diameter 
is replaced by an equivalent diameter: 

£) — 4 X till t!^CL areQ ~ (3-Ha) 

ec j wetted penmeter Ay^ 

The difference between D, and D is that D includes the wetted 

h eq eq 

perimeter of the exchanger shell. The error in assuming is 

small for most heat exchangers. 

In addition to the pressure drops due to pipe friction there are 
static pressure changes due to acceleration or deceleration of the flow 
at the heat exchanger ends and changes in specific volume. Since these 
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effects depend on pipe sizes, types of headers, etc., they are left for 
further discussion in Chapter VII. 

Correlation of Heat Transfer and Pipe Friction 

Equation (3-4b) can be expressed in terms of other dimensionless 
quantities as follows : 



y — Mu 
Re * P r 

where ^ ~ JlsL 



(3-15) 

(Nusselt No.) (3-15a) 



Pr 



— C P M 

~ K 



Re. 




(Prandtl No.) (3-15b) 



(Reynolds No.) (3-15c) 



For turbulent flow inside of 



tubes according to McAdams 



( 8 ) 



0.8 _ 0,4 

Mu = 0.02 3 R e - Pr 



(3-16) 



Equation (3-16) can be used for turbulent flow outside of tubes 
providing that the hydraulic diameter is used in lieu of d. 

From equations (3-14), (3-15) and (3-16): 



Y 



8 FV 



o^- 



(3-17) 



To summarize, the principal working equations for tubular heat 
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exchangers having longitudinal fins are: 



_L_ — _L Itj-teV 

Dk + Zl 



%x * wn 



(3-12) 



AP = r >-vi 

di 29V 



(3-13) 



Type of Heat Exchanger 

The feasibility of this cycle depends to a large extent on the size 
of the three heat exchangers; namely, regenerator, precooler and inter- 
cooler. By far the largest of these will be the regenerator. For an 
ordinary shell and circular tube regenerator analysis reveals that to 
satisfy the pressure drop and size limitations extremely small tubes 
(less than 0.2 in. diameter) placed very close together (less than 0.3 
in. center to center) are required. 

Construction difficulties, if the above dimensions were specified, 
would be of sufficient severity to render the design impractical. How- 
ever, the addition of longitudinal fins to the exterior walls of the 
tubes in this type of heat exchanger will cause a salutory effect. 
Analysis reveals that tube size and spacing may be increased, and the 
total number of tubes decreased. Although this may not be as good a 

design as could be obtained by employment of more compact surfaces such 

(4) 

as were investigated by Kays, London and Johnson the reduced complex- 
ity in design and analysis warrants its use here. 

A similar design is investigated for the precooler and intercooler, 
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I 



the difference being in dimensions. In all cases, the hotter fluid is 

outside the tubes (on the finned side). Figures 3-3 and 3-4 show the 

finned tubes and tube arrangement respectively. Although the analysis 

assumes that the fins extend the entire length, L, of the tube, boundary 

layer buildup between the fins will require the use of several shorter 

(9) 

fins staggered along the tube. This has been discussed by Gloyer 

with the conclusion that the minimum fin spacing should be not less 

than one boundary layer thickness, 

* 1 1 






-> 

-> 




Figure 3-2 



Von Raman* 6^*^ equation for turbulent boundary layer thickness, 

8 . is: 

% 



§ n 0.(183 Q 



/ \^ 5 



in. 



(3-18) 
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Figure 3-3 



Heat Exchanger Tube 





Figure 3-4 Tube Arrangement 
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For moderate Re: 




(3-19) 






_8 

7 



V, 



OV 



(3- 19a ) 



In all three heat exchangers the tubes will have six longitudinal 
BTU 

brass fins (k f ^s 60 — r ). For good heat transfer the fins 

ft - hr - °F/f t 

must have good metal-to-metal contact with the tube. The assembly 
shown in Figure 3-3 consists of two concentric tubes, the outer one hav- 
ing the six external fins. These outer tubes will be quite short (six 
to ten inches) and could be made in one machine shop operation. These 
would then have to be brazed to the inner tube in a manner to insure 
intimate contact. 

In this analysis it has been found convenient to employ areas 
associated with one tube. For the equi-angular arrangement that has 
been selected, the frontal area per tube, Cl Is equal to twice the 
area of one of the triangles formed by the centers of three adjacent 
tubes. The formulas for this and other areas are: 



A , = [ ndf .e(2«-c)\L -{f 

Ac — rrc/d it Z 

L 12 ft 



(3-20) 



(3-21) 



Af 12 Jl _ 
L 12 




(3-22) 
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( 3 - 23 ) 





* * I li- 



ft 



z 



a 



t 



S 2 Sin 60 ° fj. 2 

H + 






4 - 6 Lc 



I 

I if 




Other equations which are used in the analysis are: 



Vc. __ Q h v c 

X O c \ 



D 2 - ft 2 

ir 



__ m c y c 

0.x ox 



a 



( 3 - 24 ) 



( 3 - 25 ) 



( 3 - 26 ) 



( 3 - 27 ) 



( 3 - 28 ) 



39 






Regenerator 



/\ A, /\ 

/ \ / \ / \ 

/ \ / \ / \ 
- X / N / 

N ' \ / 

\ / ' ' 
y v 


5 






\ 


CO 

<X> 

n* 







Figure 3-5 Schematic Representation of the Regenerator 



TABLE 3-1 

PROPERTIES OF HELIUM ENTERING AND LEAVING REGENERATOR 



Location 
(Fig. 3-5 l 


rr 

o 

^3 


g (P«M) 


h (BTU/lb) 


v (ft 3 /lb) 


4 


224.9 


1000 


361.9 


1.883 


5 


874.8 


985 


1667.8 


3.677 


8 


951.3 


452.2 


1757.3 


8.413 


9 


301.5 


447.7 


951.4 


4.606 



m = 182.4 lb/sec 

At = 76.6 °F 



40 



Since this analysis is based on average values of such quantities 
as v, k, ^ , Re, Pr, Y' , and <p , they are evaluated at the 

arithmetic mean temperatures. To distinguish the hotter fluid outside 
the tubes from the cooler fluid inside the tubes the terms "hot side" 
and "cold side" will be used henceforth. 

TABLE 3-2 





AVERAGE PROPERTIES 


OF FLUIDS IN THE 


REGENERATOR 


Quantity 


Hot Side 


Cold Side 


Units 


t 


626.4 


549.8 


°F 


V 


6.510 


2.780 


ft 3 /lb 


k 


0.138 


0.132 


BTU/ft 2 - hr 




0.0759 


0.0725 


lb/hr - ft 


Pr 


0.681 


0.680 


— 


c 


1.24 


1.24 


BTU/lb - °F 



P 



Precooler 

The basic difference between the regenerator and the precooler is 
that the fluid inside the tubes for the latter is sea water rather than 
helium. The heat transmitted to the water is lost, but the lower tem- 
perature for compression results in less compressor work required, and 
because of the higher temperature differential in the regenerator, in a 
higher overall cycle efficiency. The precooler, as in the case of the 
regenerator, consists of externally finned circular tubes with single 
pass counterflow. The tubes and tube arrangement are also the same as 
in the regenerator with a slight modification of dimensions. 
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TABLE 3-3 

PROPERTIES OF FLUIDS ENTERING AND LEAVING PRECOOLER 



Location 




P (psia) 


h (BTU/lb) v 


(ft J /lb) 


i 


100 


444.4 


701.5 


3.428 


9 


301.5 


447.7 


951.4 


4.607 


in 


60 


— 


— 


.01604 


1^0 out 


90 


— 


— 


.01610 




Aj = 182.4 

At = 103 0 

TABLE 3-4 

AVERAGE PROPERTIES OF FLUIDS 


lb/sec 

F 

IN THE P RE COOLER 


* 


Quantity 


Helium 


h 2 o 


Units 




t 


200.8 


75 


°F 




V 


4.018 


0.01607 


ft 3 /lb 




h 


0.0558 


2.22 


lb/hr - ft 




c 

P 


1.24 


1.0 


BTU/lb - °F 
? 




k 


0.10 


0.346 


BTU/ ft" - hr 


- °F / ft 


Pr 


0.692 


6.42 


— 
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Intercooler 




Figure 3-7 Schematic Representation of the Intercooler 



TABLE 3-5 

PROPERTIES OF FLUIDS ENTERING AND LEAVING INTERCOOLER 



Location 


t_cn 


P-Cpsial 


h (BTU/lb) 


v (ft 3 /lb) 


2 


216.9 


666.7 


848.8 


2.768 


3 


100 


663.4 


703.7 


2.309 


R 2 0 in 


60 


— 


— 


0.01604 


11^0 out 


90 


— 


— 


0.01610 




“h 


= 182.4 


lb/sec 






At 


= 75.4 


°F 








TABLE 3- 


6 






AVERAGE PROPERTIES OF FLUIDS IN THE INTERCOOLER 


Quantity 


Helium 




h 2 o 


Units 


t 


158.5 




75 


°F 


V 


2.538 




0.01607 


ft 3 /lb 


h 


0.0532 




2.22 


lb/hr - ft 


c 


1.24 




1.0 


BTU/lb - °F 


p 










k 


0.096 




0.346 


BTU/ft 2 - hr - 


Pr 


0.687 




6.42 


— 



°F/ft 
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Selection of Tube Dimens i o ns 



The dimensions that appear in Table 3-7 are the result of many 
trial and error calculations. There are, no doubt, other combinations 
which will produce equally good results. It is intended that the ones 
used here represent a reasonable design for this cycle. 

Calculations 

Values in Tables 3-8 and 3-9 were obtained from substitution of 
known and assumed quantities into equations outlined earlier in this 
chapter. It can be seen that there occur several relations and vari- 
ables. However, there is only one independent variable and it is 
convenient to choose as this variable. Figures 3-8, 3-9, and 3-10 
show the variation of AP, L, D and total volume with for the three 
units . 

The final selection of V^depends on how large a pressure drop 
and/or size can be tolerated. The choices of V, and the corresponding 

n 

heat exchanger dimensions listed in Table 3-11 are based on the values 
assumed for AP/P. in Chapter II. A knowledge of end effects was 
necessary to make the selection since in some cases the AP due to 
the end effects is larger than AP due to pipe friction. This is 
covered in more detail in Chapter VII. 
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TABLE 3-7 
TUBE DIMENSIONS 



' 


Regenerator 


Intercooler 


Precooler 


Units 


d. 

1 


0.30 


0.40 


0.40 


in. 


d 

o 


0.35 


0.50 


0.50 


ii 


d f 


0.40 


0.60 


0.60 


it 


A 


0.20 


0.20 


0.20 


M 


c 


0.025 


0.05 


0.05 


M 


s 


0.80 


1. 00 


1. 00 


11 







TABLE 3-8 








HEAT EXCHANGER 


AREAS AND ASSOCIATED 


QUANTITIES 






Regenerator 


Intercooler 


Precooler 


Units 


Vl 


0.2922 


0.3321 


0.3321 


ft 2 /ft 


A c/L 


0.07854 


0.1047 


0.1047 


11 


A f/L 


0.20 


0.20 


0.20 


II 


Cl h 


0.00277 


0.00363 


0.00363 


ft 2 




0.000491 


0.000873 


0.000873 


11 


Qt 


0.003849 


0.006014 


0.006014 


II 


D h 


0.0379 


0.0437 


0.0437 


ft 


£ 


0.1774 


0.2405 


0.2405 


— 


n 


0.2688 


0.3152 


0.3152 


— 


w 


1.0 


3.897 


6.717 


— 
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TABLE 3-9 



LIST OF CONSTANTS FOR THE REGENERATOR, FRECOOLER AND INTERCOOLER 



Quantity * 


Regenerator 


Precooler 


Intercooler 


Reference 

Equation 


Vc/V h 


2.407 


0.1338 


0.127 


(3-26) 


Re h/v h 


276.1 


700 


1165 


(3-15c) 


Re At 

c/V 

c 


446 .5 


3360 


3360 


(3-15c) 


hh/ n v h 


685.7 


1110 


1759 


(3 -4b) 


zr 

o 

o 

< 

o 


1605.7 


224,000 


224,000 


(3 -4b) 




0.1575 


0.1558 


0.1566 


(3-17) 


w/ c 


0.1575 


0.04098 


0.04098 


(3-17) 


ra ^//ii h 


0.0667 


0.0471 


0.0471 


(3-10a) 


At>h/L 4 v h 


4.365 x 10" 4 


6.14 x 10“ 4 


9.73 x 10' 4 


(3-13) 


APc/L ?c V c 


15.53 x 10' 4 


2010 x 10‘ 4 


2010 x 10' 4 


(3-13) 


L *K 


0.0805 


0.02137 


0.01694 


(3-12) 


d 2 /n 


0.00490 


0.00765 


0.00765 


(3-27) 


NV, 

h 


4.288 x 10 5 


2.02 x 10 5 


1.2753 x 10 5 


(3-28) 



* alL units are as indicated in the list of 
symbols at the , beg inning of the Chapter. 
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■’igure 3-9 

___J ' L 



Variation tof frecoolcr Characteristics with Average 
;telUua~V«lo«jLty 




Total p olurae (ft'3 



TABLE 3-10 



DATA FOR SELECTED HEAT EXCHANGERS 





Regenerator 


Precooler 


Intercooler 


Units 


v h 


41.5 


75 


75 


ft/sec 


V 


100 


10.4 


9.5 


M 


c 


L 


35.7 


9.9 


10.2 


ft 


D 


7.19 


4.52 


3.61 


it 


N 


10,330 


2690 


1700 


Tubes 


Volume 


1450 


161 


104 


ft 3 


* p h 


0.76 


0.715 


1.05 


psi 


AP 


12.03 


4.67 


4.17 


t ? 


c 


a 


10.6 


21.4 


24.5 


in. 
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CHAPTER IV 



COMPRESSORS 



Symbols used in this chapter: 

Symbol Description 

A annulus area 

C force coefficient 

c chord 

c specific heat at constant pressure 

. Last stage isentropic enthalpy change 
ra ° First stage M " M 

d diameter 

F blade force 

f reheat factor 

g conversion factor, 32.17 



H 



total enthalpy 
static enthalpy 



or 



blade height 

J conversion factor, 778 

last stage actual enthalpy change 
L ratio first „ 

m mass flow 

N rotational speed 

n no. of stages 

P pressure 

r radius 



Units 

2 

in. 



in. 

BTU/lb - °F 



in . 



lb 



f 



lb^ sec 1 
BTU/lb 
BTU/lb 
in . 

ft - lb 
BTU 



lb/sec 

rpm 

psia 

in. 
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Symbol 



T 

U 

V 

v 

W 

oc 

/3 

A 

€ 

P 

<P 

X 

co 



Description 
entropy 
blade spacing 
absolute temperature 
peripheral velocity 
absolute velocity 
specific volume 
relative velocity 
absolute flow angle 
relative flow angle 
increment 

drag force/lift force 

efficiency 

density 

defined by equation (4-8) 
defined by equation (4-7) 
degree of reaction 
rotational speed 



Subscript Description 
C compressor 

c center, mean 

D drag 

DP drag profile 

e exit 

i inlet, initial 



Units 

BTU/lb - °F 
in. 

°R 

ft /sec 
f t/sec 
ft 3 /lb 
ft /sec 
degrees 
degrees 



lbs/ft 3 



rad/sec 
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* 







Subscript Description 



L 



lift 



ra 



axial 



n 



last stage 



R 



rotor, relative 



S 



stator 



s 



stage, isentropic 



th 



theoretical 



u 



tangential 



x 



arbitrary stage 



General 



Before proceeding with a design study of the two compressors, the 
type of compressors must be selected. From Table 2-2 the energy impart- 
ed to the gas in the two compressors is 147.3 and 158.2 BTU/lb respect- 
ively. Analysis reveals that each compressor, if centrifugal, would 
consist of two or more stages in order to prevent excessive peripheral 
speeds. Rather than employ several centrifugal stages where efficien- 
cies are usually low due to losses occurring between stages, two axial 
flow compressors have been selected. 

The flow channel associated with an axial flow compressor is 
illustrated in Figure 4-1 and the corresponding h-s diagram in Figure 
4-2. 



53 



Inducer 




u 

o 

w 

cn 

<D 

U 

B 

o 

o 

5 

o 

r~< 

tx* 



X! 

< 

C 

C3 

M-J 

O 

c 

o 

*r 4 



CO 

OJ 

c: 

Q) 

CO 

0) 

v-i 

a 

<u 

od 

o 

•r-< 

-u 

§ 

0) 

JC 

o 

CO 



<r 

<u 

3 

to 



54 




Figure 4-2 Typical Enthalpy -Entropy Diagram for an Axial 

Flow Compressor 
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Prom Figure 4-2, where it la aseSned that the flow in the inducer 
and diffuser sections is steady and adiabatic, the following equations 
are obtained if geopotential energy differences are ignored: 



H = 


V 2 

2 9 J + h 


* constant 


(4-1) 


^5 

r • 

III 






(4-2) 














* 7 e — 


(a h s ) nj e 




(4-3) 







(4 -4a) 



Compressor Stage Performance 

Before the compressor design carxjbe carried out the type of stage 
must be selected. In this design only one type of stage was considered; 
namely, a non-symmetric stage where the absolute velocity entering the 
rotor is axial. Figures 4-3 and 4-4 show this type of stage and a 
corresponding h~s diagram. 
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Stator 




I 





Figure 4-3 Schematic Representation of an Axial Flow Compressor 

Stage and Corresponding Velocity Diagram 
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Figure 4-4 Typical h-s diagram for an axial flow compressor stage 



For the purpose of developing convenient working equations for 
calculating stage performance the velocity diagram of Figure 4-3 is 
changed to a dimensionless form by dividing all vectors by U to obtain 
Figure 4-5. It is assumed here that • 



58 




Figure 4-5 Dimensionless velocity diagram for axial flow 

compressor stage. 



In Figure 4-5 several new terms have been introduced. They are: 






R 



Drag Force 
Lift Force 



for rotor blade 



(4-5) 



Drag Force 

$ Lift Force for stator blade 



(4-5a) 






W u , + Wu a 

2U 



T 

2 



(4-6) 
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(4-7) 



<p = = ~ 

r U U 






V < 2 ~ 

u 



Vy^ 2 - Wu, 

u 



(4-8) 



or 



_ V, 



42 



u 



for oC = O 



(4-8a ) 



Expressions for stage efficiency are based on airfoil theory 

assuming incompressible flow and will not be derived here. Referring 

to Figure 4-5; if the product 6 6 c * s Ignored with respect to £ 

Kb * 



or ^s* 



P*-/? 



_ U ±C 

vs 



+ X - £r 4> 



<f> + £ s (i-x) <P + £r * 



(4-9) 



uV 



fe vs 



(f* = e S=0) <*-9a) 



% 



fi-f? 



4 > 



/ X 6s 4* X — ^ 

(p + e s (i-$ <p + <f>x 



(4-10) 



where £ — 



^4 



(4-5) 
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+ o.oi a C Z L 



(4-11) 



Q = C DP + ao! °i 



c L = 

C o 



_ 2t(«, 



Wo /u 



(4-12) 



C, = 



2r(£X 






/u 



(4-12a) 



Equation (4-11) was obtained from Howell^^ and is based on 
cascade tests and experimental single-stage and multiple-stage compres- 
sor data. These tests also show that the profile drag coefficient, 

Cpp, is primarily a function of the ratio s/c. This relation is plotted 
in Figure 4-6. (Howell) 

Figures 4-7 and 4-8 contain plots of additional test results^^, 
showing a relation between fluid deflection, exit flow angle and ratio s/c. 

In addition to the stage efficiency, calculation of energy input 
per stage is necessary for determining the number of stages required. 

This can be expressed as: 




(4-13) 
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<AHj= (&H) * % 



(4-14) 



U r ns 
9 J 



For a Vy 



^s)s/ 



Uti), 

9 J 



(4 -14a) 



The degree of reaction is defined as the ratio of actual enthalpy 
rise in the rotor to actual enthalpy rise in the stage. 



Degree of reaction 



— hz zhi — h-Lihi 

h 3-h, h 3 - h, ~ (ah) 



(4-15) 



3,1 



If the relative flow in the rotor blade is steady and adiabatic: 



Hr — h + 



u 



W 



2 9 J - z. g J 



=3 constant (4-16) 



For U- = U 0 and W * W : 
12 



k 2 -k, — 



W | - W2 

2 a J 



2 l 

W U| - 

2 0 J 



(4-17) 
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Proa (4-8) , (4-13), (4-15) and (4-17): 



Degree of reaction = 



W-t, + W a 2 
2 U 



X 



(4-6) 



Determining Number of Stages 

In a multistage compressor, the stage performance varies from 
stage to stage in some continuous manner. By evaluating the perform- 
ance of the first and last stages where the stream properties are known 
and assuming some reasonable relation between these and the other 
stages' it is possible to express U, X and for any stage in terms 

of the two known stages. 

The relations assumed in the analysis are: 




x-j 



(4-18) 




(4-19) 



x-i 




(4-20) 



From f4-l4a) it follows that: 




(4-21) 
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(4-22) 



X ( Ah s)x — (/ + f)(AHs) /n — (Ahs), I.D^D 
X=l L * = * 



~ *“] 



where 



P) = (Ahsln / Un\ TTn Hn 

(Ahs), 'W L i Hi 



or 



0 + ^(4 jj§kn == n = | + D + 



Qzi \ 

D ‘ ‘ 



n-i 



D fR-| 



(4-22a) 



(4-22b ) 



where f) * number of stages required if for 
all stages (Ah g ) x = (dh^j 



Similarly: 




(4-23) 








(AH)n 

(Ah), 



n !i = I + L + L 

„here I = ^2 

(Ah), 




(4-24) 



(4-25) 



(4-25a) 



and = Number uf stages required if for 

all stages (Ah)^ = (Ah)^ 
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Since there must be an integral number of stages the solution of 
equation (4-22b) for (&h ), generally will not be the same as the 

S 1. 

value obtained from the first stage calculations. Thus a correction to 
the peripheral speed is usually necessary. Unless this correction is 
large, the effect on stage efficiency will be small. The total energy 
input can then be calculated with equation (4-25). 

The overall compressor efficiency is defined as: 



1 t 



(a Hs) i ; e 

(AH)- S e 



(4-26) 



where 






(l-6a) 



and 





(4 -4a) 



The reheat factor, f t appears in equation (4-22b)« Equations (1-8) 
and (1-9) give the reheat factors for an infinite number of stages for 
both expansion and compression processes. Referring to Figure 4-9 the 
reheat factor for H stages of compression is obtained from: 

n 

(Ah 

J 

Hf'-hi. 




(4-27) 
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Figure 4-9 



For an infinite number of stages (4-27) becomes 



l + f 



co 



oc 

E 

X~t 

h t ' - hi 



(4-27a) 



Since dh = Tds + vdp for either a reversible or irreversible 
process and dh = Tds when P = constant: 



Ah) — I To/s = area under a constant pressure 
k /p__ r J 

on * line on T-*S diagram 
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( Ah s)x 



area between the stage inlet and 
exit constant pressure lines on the 
T-S diagram extending to a reference 
absolute temperature . 



r dh f 

The difference between and ^ n s^x is the small 

Xrl Xnl 

triangular areas shown in Figure 4-9, Assuming the constant pressure 

lines to be parallel and ( Ah ) constant for all stages, the area of 

s x 

2 

i-f'-f consists of n small triangles. 

Equations (4-27) and(4-27a) can be combined as follows: 



h+'-hi ~ (^s)x 

T Xf=J 

hf'-hi - £ dh s 



or *f 

K 



area i-f-f 1 less n small triangles 
area i-f-f 



f 




n-i 



n 2 



n 



f Od (4-23) 

'oo n 



Compressor No. 1 

Having previously decided on the type of blading it is now neces- 
sary to determine an optimum combination of the many variables involved. 
First of all, a high stage efficiency is desired. To calculate this, 
practical values for (s/c)„, (h/c). and (h/c) must be selected from 

Kb K 
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experience. The quantities, C Dp , (s/c) g , and A/3 , can be obtained 
from test data of Figures 4-6, 4-7 and 4-8. Then for given values of 
/S 2 the remaining quantities in the dimensionless velocity diagram, 
Figure 4-5, and the stage efficiency can be calculated. 

The effect of a bad choice of (h/c)_ and (h/c) is not great, 

S R 

since these quantities change the stage efficiency slightly. The 
values assumed for a preliminary investigation were 3.0 and 2.0 respect- 
ively. The choice of (s/c) p is more critical. From a consideration of 
dimensions, especially blade heights, (s/c) of less than 1.0 proves to 
be impractical. Therefore, stage efficiencies were calculated for 
(s/c) R of 1.0 and 1.5. Figure 4rl0 shows the plot of g vs. ^3 ^ f° r 
these two ratios. It can be seen that the maximum efficiency in both 
cases will occur at some smaller angle than has been plotted. At the 
smaller angles the blade heights become smaller than practical design 
permits. In addition, consideration must be given to Compressor No. 2 
and Turbine No. 1 since their rpm will be the same as Compressor No. 1. 
For the same blading in Compressor No. 2 the smaller specific volume 
will result in even smaller blade heights. Too high a rpm will result 
in a large through flow velocity in Turbine No. 1 and consequently a 
low turbine efficiency. 

In conventional designs using air as the working fluid, the maxi- 
mum peripheral speed is usually determined from Mach effects because of 
the relatively low velocity of sound in air. The velocity of sound in 
helium is approximately three times as great so that allowable stress 
becomes the controlling factor. The materials that are presently being 
used for compressor rotors and blades will permit a U max of about 
1200 ft/sec. 
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Figure 4-11 shows plots of N vs . /Q for last stage blade heights 
of 1.0 and 1.4 inches. For calculation of the last stage blade heights 
it is assumed that the same blading is employed as in the first stage. 
From the plot it can be deduced that for given values of and 

the compressor having the smallest value of (s/c) has the lowest rpm. 
The compressors in this cycle require high rpm if small blade heights 
and large diameters are to be avoided. High turbine rpm results in a 
blade height to mean diameter ratio which is too large, so some compro- 
mise has to be made. Preliminary calculations indicate that 7000 rpra 
is satisfactory for Turbine No. 1 and the two compressors. Another 
factor "to be considered is the number of stages required for different 
values of (s/c)^, 

TABLE 4-1 



EFFECT OF (s/c ) 0 ON COMPRESSOR BLADING 
K 



(S/ C) R 


1.0 


1.5 


Units 


N 


7000 


7000 


rpm 


Si 


62.2 


64.4 


degrees 


h 


1.4 


1.4 k 


inches 


n 








> 


.882 


.873 


— 


n 


7 or 8 


10 


stages 
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From the comparison in Table 4-1 the blading having a space to 
chord ratio of 1.0 appears to be the most favorable and is used in 
subsequent analysis. 

The results of preliminary calculations indicate that an (h/c) 

R 

or 3 results in a very small chord so it was changed to 2.0. 

In order to make all rotor disks of the same dimensions, the blade 
root diameter will be constant. The first and last stage data that 
result are listed in Table 4-2. 

TABLE 4-2 

DATA FOR FIRST AND LAST STAGES OF COMPRESSOR NO. 1 





First Stage 


Last Stage 


Units 




0 


0 


degrees 




1.0 


1.0 


— 


ftfl* 


2.0 


1.63* 


— 


( h /c) S 


2.0 


1.63* 


— 


ft 


62.0 


62.0 


degrees 


( 3 , 


70.6 


70.6 




<P 


0.3522 


0.3522 


— 


r 


0.3376 


0.3376 


— 


X 


0.8312 


0.8312 


— 




43.8 


43.8 


degrees 




0.48 


0.48 


— 


% 


0.8812 


0.8755 


— 




0.0482 


0.0513 


— 




0.0438 


0.0440 


hn 



* h/c for last stage = h/c for first stage x — 
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Since there must be an integral number of stages the calculation 

of work imput per stage, dimensions, and velocities requires a series 

of approximations. Table 4-3 contain^ firs t approximations of U, V^, 

A , h and d^ f for a fixed compressor speed of 7000 rpm. Before 

proceeding with the analysis, some additional information is necessary; 

namely, the inlet and exit pipe velocities and the inducer and diffuser 

efficiencies. The velocities were selected after taking into consider- 

' * 

ation the size of pipe required and the efficiencies are typical values 
for conventional designs. These are indicated in Table 4-4 along with 
calculated results. 

The following equations, supplementary to those listed previously, 
are also applicable: 



A ~ — - x M4- — t rhd c (4-29) 

V" 

N — rpm (4-30) 

Tide 

TABLE 4-3 

FIRST APPROXIMATION DATA FOR FIRST AND LAST STAGES OF COMPRESSOR NO. 1 





First Stage 


Last Stage 


Units 


N 


7000 


7000 


rpm 


u* 


1200 


1190 


ft/sec 


V 


423 


420 


f t/sec 


A * 


213 


172 


. 2 

in 


h* 


1 .72 


1.40 


in 


d * 


39.32 


39.00 


in 



c 



* First approximation 
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TABLE 4-4 



ADDITIONAL DATA FOR COMPRESSOR NO. I 



v i 


400 


ft/sec 


V 

e 


400 


tl 


nt 


0.75 


— 


1e 


0.85 


— 




124.1 


BTU/lb 


m,; 


124.3 


It 


(& h 3 ) ( * 


17.14 


tt 


o’* 


7.34 


— 


i+f 


1.0117 


— 


n 


7 


stages 


D 


0.9773 


— 


L 


0.9837 


— 


n‘ 


6.920 


— 


n" 


6.941 


— 




18.17 


BTU/lb 


(An), 


20.60 


ti 




143.05 


tt 


% 


0.868 


— 


% 


0.867 


— 



* First approximation 
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We observe in Table 4-4 that the actual required for seven 

stages differs from the first approximation. Therefore, some changes 
must be made in U and some of the dimensions if N is to remain at 7000 
rpm. These changes can be made by use of the following equations: 



I 




(4-31) 



d 



* u 
c u 



(4-32) 



X = v * -Jp (‘-w 

A * I/ 

A, = A, -jj- <*-*> 



Corrected results for the first and last stages are indicated in 
Table 4-5. 



76 



TABLE 4-5 



FINAL DATA FOR FIRST AND LAST STAGES OF COMPRESSOR NO* 1 





First Stage 


Last Stage 


Units 


u 


1235 


1225 


ft/sec 


N 


7000 


7000 


rpm 


V 1 


436 


432 


ft/sec 


V 2 


604 


599 


II 


W 1 


1310 


1300 


If 


W 2 


928 


920 


It 


A 


207 


167 


. 2 
m 


h 


1.63 


1.33 


in 


d i 


38.87 


38.87 


ii 


d 


40.50 


40.20 


n 


c 


d 


42.13 


41.53 




o 


c 


0.82 


0.82 


n 


S R 


0.82 


0.82 


ii 


S s 


0.39 


0.39 


ii 



Compressor No. 2 

The same blading is used in Compressor No* 2 as was used in Com- 
pressor No. 1 and the rpm is fixed at 7000. Figure 4-12 shows the 
variation of h^ and n as a function of U. A mean peripheral speed of 
1050 ft/sec was selected in order to obtain a satisfactory last stage 
blade height. The root diameter is constant as in Compressor No. 1. 

Table 4-2 is applicable with the following exceptions for the 
last stage: 
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(f>/c) R = ( h /c)s = 166 

rj $ — O.Q-76Z 



Calculated results are listed in Tables 4-6, 4-7 and 4-8. 



TABLE 4-6 

FIRST APPROXIMATION DATA FOR FIRST AND LAST STAGES OF COMPRESSOR NO. 2 





First StaRe 


Last StaRe 


Units 


u* 


1050 


1037 


ft/sec 


N 


7000 


7000 


rpm 


V 


3 70 


365 


ft/sec 


A * 


164 


134 


. 2 
m 


h* 


1.52 


1.26 


in. 


d * 


34.38 


33.90 


n 



c 



* First Approximations 
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TABLE 4-7 



ADDITIONAL DATA FOR COMPRESSOR NO. 2 



v i 


300 


ft/sec 


V 


450 


it 


e 






Hi 


0.75 


— 


Oe 


0.75 


— 




126.6 


BTU/lb 


(AHtfn 


128.83 


»t 


fast* 


13.12 


m 


n'* 


9.97 




1+f 


1.0157 


— 


n 


10 


stages 


D 


0.9793 


— 


L 


0.9849 


— 


n' 


9.9005 


— 


n" 


9.9189 


— 


(4 


13.20 


BTU/lb 


(Ah), 


14.98 


n 


(* Hl.ii 


148.58 


ii 


'll,. 


.867 


— 


1 ' 


.840 


— 
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TABLE 4-8 



FINAL DATA FOR FIRST AND LAST STAGES OF COMPRESSOR NO. 2 



First Stage 



Last Stage Units 



U 

N 

V 1 

V 2 

W 1 

W 2 

A 

h 



d. 

i 



c 

s. 



1050 

7000 

370 

513 

1115 

788 

164 

1.52 

32.86 

34.38 

36.90 

0.76 

0.76 

0.37 



1037 


ft/sec 


7000 


rpm 


365 


ft/sec 


507 


n 


1100 


IT 


778 


II 


134 


in 2 


1.26 


in* 


32.64 


it 


33.90 


1 1 


35.16 


it 


0.76 


n 


0.76 


i« 


0.37 


i? 



r 
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CHAPTER V 



TURBINES 



Symbols used in this chapter: 



Symbol 


Description 


Units 


d 


diameter 


in. 


f 


reheat factor 


— 


8 


conversion factor, 32.17 


lb ft 

m 






lb f sec 


H 


total enthalpy 


BTU/lb 


h 


blade height 


in. 


or 


enthalpy 


BTU/lb 


J 


conversion factor, 778 


ft - lb 
BTU 


K 


coefficient defined in equation (5-10b) 


— 


m 


mass flow 


ibs/sec 


N 


rotational speed 


rpm 


n 


• 

number of stages 


— 


P 


absolute pressure 


psia 


r 


radius 


in. 


T 


absolute temperature 


°R 


U 


peripheral speed 


ft/sec 


V 


absolute flow velocity 


ft/sec 


V 


specific volume 


ft 3 /lb 


W 


relative flow velocity 


ft/sec 


oC 


absolute flow angle 


degrees 


P 


relative flow angle 


degrees 
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Symbol 


Description 


Units 


A 


increment 


— 


n 


efficiency 


— 


X 


defined by equation (5-16) 


— 


V 


defined by equation (5-15) 


— 


<$> 


coefficient defined by equation (5-7) 


— 


X 


defined by equation (5-20) 


— 


X 


degree of reaction 


— 


r 


coefficient defined by equations (5-12) 
and (5-13) 


— 


Od 


rotational speed 


rad/sec 


Subscript 


Description 




B 


blade 




c 


center, mean 




e 


exit 




i 


inlet, initial 




m 


axial 




n 


last stage 




o 


outer 




R 


relative 




s 


isentropic 




T 


turbine 




u 


tangential 




General 




This 


chapter consists of a design study of the 


two turbines in the 


cycle. As 


in Chapter IV for the compressors, a one* 


-dimensional flow 
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analysis is performed to determine turbine diameters, blade heights, 
number of stages and blade shapes at a mean streamline. The flow is 
assumed to be steady and adiabatic in both moving and stationary 
channels so that the total enthalpy, H, in the stationary channels and 
the total relative enthalpy, in the moving channels are constant 

along the streamlines. Analysis reveals that if a centripetal turbine 
is employed an excessive peripheral speed will be required to produce 
the necessary work output. Therefore, the entire analysis that follows 
is based upon an axial flow turbine. 



Stator 




Figure 5-1 Schematic Representation of an Axial Flow Turbine 



34 



Figure 5-1 is a sketch of an axial flow turbine showing stream 
locations 1, 1, n and e. The inlet and exit sections serve to deceler- 
ate and accelerate the flow respectively, since, as we shall see later, 
the pipe velocities for this cycle are necessarily larger than the 
through flow velocity. Referring to the enthalpy - entropy diagram in 
Figure 5-2, the following relations can be written: 




where 




(5-la) 



and 




(5 -lb) 



Hi — H, and He — Hn 



(5-2) 



where 




(5 -2a ) 



AH e 



H,-H n 



(5-3) 




(5-4) 
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Figure 5-2 Typical Enthalpy-Entropy Diagram for an Axial Flow Turbine 



Stage Efficiency 

Figures 5-4 and 5-5 are enthalpy - entropy diagrams illustrating 
absolute flow quantities for the stage and relative flow quantities for 
the rotor blade respectively. 

From Figure 5-4: 



(Alls) 


= h,» - h 3 > 


(5-5) 


Kf 




(5-6) 


/ 

(Ah,) 


— ^ 2 * h 3?^ 


(5-7) 


Assuming zero reheat (f « 0): 




ii 


= 


(5 -6a) 


4> 2 = 


Ah),;* 


(5-8) 




(Ahs)” . „ . , 

-i s isentropic degree 

of reaction 


(5-9) 


III 

10 

c- 


Ah B — Ahe v > 

(Ah s V (Ahs)” 


(5-10) 
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Stator 



/////// / 



/ 




2 

« • 




3 











F\ ' \ \ \ \ \ \ \ 



Rotor 




(a) 



(b) 



Figure 5-3 Velocity Diagram for an Axial Flow Turbine Stage 







Figure 5-4 Typical Enthalpy-Entropy Diagram for an Axial 

Flow Turbine Stage 



88 



In equation (5*10) the kinetic energy of the fluid at the exit of 
the stage is not considered as useful energy. Actually, most of this 
energy is available for the succeeding stage and another stage efficien- 
cy is defined: 




(5-i0a) 



where ( A h s f == + K 



(5-10b) 



and K * percentage of the kinetic energy 
leaving a stage that is available 
for useful work in the succeeding 
stage 



89 




Figure 5-5 Typical Enthalpy - Entropy Diagram for Relative Flow 

in a Turbine Rotor Blade 



From Figures 5-4 and 5-5: 



H Rj - constant 




(5-11) 



(5-lla) 



(5-12) 
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(5-13) 



% 



W3 

W 3 ’” 



Wg \y 2 

2 9 J ~ ^ 






(Ah s ) 



+ 



2g J 



(5-14) 



For zero reheat (f = 0): 



w , 2 _ 



U , 2 



= Vi 



2 q T (Ah s )’x ’ ^ V^w/t U3 






u, 



(5 -14a) 



It is advantageous for subsequent analysis to introduce the 
following definitions: 



y 



u; 



( Ah s)’ ( 2 3 J ) 



(5-15) 






.. a 



U; 



(5-16) 



V, 



^3 



V, 



m 2 



(5-17) 



From the trigonometry of Figure 5-3, the definitions above and 
equation (5-14a), several new equations can be derived: 
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(5-14b) 



X’ f 4>* %*( I-/') - 
Z (j> ii 2 ]/V r l-X’ Sin oC^ + +A 2 -lj 

^ — X - ^ 2 . X (5-18) 

U 2 A 3/ A 

*J 5 = £ V (<pvTlC Sin*C t + X" >? v) (5 ' l9) 

X* = X ; Cos* oC t )(|-X*) - 

2 $ ^ Vl-X? S/n <* 2 + 7 ( 5 ' 2 °) 




to n /&, 



ton oc i 



V 

^l/Px 5 



(5-21) 



The work output per stage can also be expressed in the familiar 



form: 



* h B= (* t v„- R,\g 

if R, = R 2 : 

All 0 = Ifj ( V ^2 ” K 3 ) 
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(5-22b) 



If = 0 ■ 



Ah 



B 



UVu, 

9 J 



When all stages are similar and produce equal work: 



AHb = n Ah B 



(5~22c) 



From the definition of reheat factor and equation (5-10b): 



A u \ = _ 

l+f “ 



h Ah, 






1* 



2g J 



+f 



(5-23) 



where 



f - f. 



n-i 

n 



(4-28) 



Two additional relations necessary for turbine design are: 



7 20 U 



I _ tcu 

Gc 7T N 



tn. 



(5-24) 



A = TT d c h — ^ XIH- Ln z 



\4 



(5-25) 



m 



Selection of Blading for Turbine No. 1 

The type of blading to be selected must permit a good efficiency 
as well as reasonable dimensions and number of stages for the condi- 
tions imposed. The use of low efficiency Curtis impulse stages would 
not be satisfactory despite the fewer stages required. Generally, it 
can be expected that reaction stages will have the highest efficiency 
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and require the greatest number of stages. In order to obtain a high 
efficiency and prevent too many stages a compromise between impulse and 
reaction stages may be warranted. Two general types of stages were 
investigated; namely, (1) symmetric or 50 per cent reaction and (2) 
modified Rateau. The modified Rateau stage differs from the pure 
impulse stage in that a small amount of reaction is introduced by 
specifying an axial exit velocity. 

Comparisons of the two types have been based on conditions exist- 
ing at Turbine No. 1, the rotational speed of which is 7000 rpm. The 
maximum peripheral speed, U, is limited by stress conditions. At 
temperatures of 1350°F a peripheral speed of about 800 ft/sec seems 
possible with materials currently available. Losses in the stator and 
rotor, expressed by the definitions of (p , and are est *~ 

mates based on conventional gas turbine performance tests. For a 
single stage the differences between and can be considered negli- 
gible. In addition, the tTTrough-f low velocity, V^, is assumed to be 
constant . 

For each value of there is an efficiency curve which can be 

plotted as a function of — - . The choice of V depends to a large 

v 2 m 

extent on the ratio of blade height to mean diameter. From Table 2-2: 

AH g = 305.9 BTU/lb 
zzz. 182.4 Ib/sec 

V, = 5.019 ft 3 /l b 
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From equations (5-24) and (5-25): 



cL = Z 7.(o 5 
s 

= I 51 7.6 {. t .z 1 .0 

Sec 

V » 300 ft/sec results in a first stage —r~ ratio of 0.18. Velo- 
m d 

c 

cities less than 300 ft/sec result in higher efficiencies but the blade 
heights become too large, especially at the last stage. 




(a) Modified Rateau Stage 




Figure 5*6 
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Using the quantities assumed below, the number of stages, n, and 
stage efficiency, ^ , were calculated for the two types of stages 
indicated in Figure 5-6. 

N = 7000 rpm Cp = 0.95 

U = 800 ft/sec y/ = 0.95 



/U = 1.0 0.92 

A - 1.0 K = 0.90 



Plots of and n vs. • are shown in Figures 5-7 and 5-S. 



It is apparent that there is little difference in maximum efficiency 
and that two less stages are required for the modified Rateau stages. 
For this reason the modified Rateau stage was selected. At the condi- 
tion of maximum efficiency: 




Ah 6 - 35.3 BTU/lb 

■7J = 0.854 

X * 0.333 



n « 8.7 theoretical stages 

An adjustment must be made in the work per stage to increase the 
number of stages from 8.7 to 9. For a nine stage turbine having equal 
work per stage: 






_ 34.0 gTU/lb 
9 



Since 




UVut _ u2 /VO 
e cr “ g 7 \ u/ 



(5-22b) 
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(Ij 2 ) — Ui(Vu a /Ua) 

\Jz)actu«l — Tw" 7 n 

(Vc, 2/ / Uj actual 



X 



(Aba) actual 

A h e 



To use the same blading, — = constant. 

U 2 

* l\ i 2 \ — ( i z „ C^^e)4ciua.l 

• • iyiLt ««i — u * x ^rfr 



i 

= (m) x 800 = 785 /t/sft 



(5-26) 



Recalculation of velocities and dimensions assuming that the 
blade root diameter, d^, is constant along the turbine, produces the 
results contained in Table 5-1. 

In calculating the turbine efficiency, inlet and exit effects 

and reheat factor must be known. The reheat factor can be calculated 
by the relation: 



{ = fco <*-*•> 

where foo is obtained £rom equation (1-9). 

To obtain (A Hs\ e ’ the values of v e > ^ and ^ G rau3t 
be known. The pipe velocities that are assumed here will be covered 
i»i more detail in Chapter VII. 



V. 

1 


■ 600 


ft/sec 


V 

e 


= 500 


ft/sec 




* 0.85 




0 e 


= 0.75 
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Using these data and equations (5-23), (5-1) and (5-4): 



1+f = 


1.0079 


(aHs)^- 


356.9 BTU/lb 


(kHc^ e = 


356.3 BTU/lb 


1 T = 


0.858 



Recalculation of the exit pressure and specific volume is not 
necessary since the calculated efficiency is only slightly higher than 
the assumed value. A summary of the first and last stage data for 
Turbine No. 1 is listed in Table 5-1. 
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TABLE 5-1 





FIRST AND LAST 


STAGE DATA FOR TURBINE NO. 


1 




First Stage 


Last Stage 


Units 


V 2 


1121 


1162 


f t/sec 


V 3 


294 


305 


ft 


*2 


74.78 


74.73 


degrees 


<*3 


0 


0 


ii 


W 2 


418 


433 


ft/sec 


W 3 


847 


878 


if 


^2 


45.29 


45.29 


degrees 


@3 


-67.98 


-67.98 


it 


U 


785 


313 


ft/sec 


d. 

i 


21.89 


21.89 


in. 


d 

c 


27.15 


28.10 


u 


d 

o 


32.41 


34.31 


tt 


h 


5.26 


6.21 


tt 


h/d 

c 


.194 


.221 


— 



Additional Turbine Data : 

N » 7000 rpm 

7$ - 0 . 854 

% = 0.858 

n = 9 stages 
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Turbine No. 2 



Unlike Turbine No. 1, this turbine is not restricted to a specific 
design rpm, thus the speed of this turbine is an additional variable. 
Since the conditions for this turbine are not very much different from 
those of Turbine No. 1, a similar type blading is used for simplicity. 
Thus, U, V^, ^ , <p , y 2> and K remain the same. From 

Table 2-2: 



AH b - 193.7 BTU/lb 

» 182.4 lbs/sec 

SJ X « 6.692 ft 3 /lb 



Since 



720 U 
i r Nl 



i.n 



(5-24) 



and 



A =7Td c h 




x 144 



un 



(5-25) 



h — n'l.Vj N * 144- . 0.001017 N Inches 

' v m , U, * no 

For N » 5500 rpm: 

a 5.59 inches 
^ c ^ s 33.33 inches 

h l/d a 0.168 

c. 



The theoretical number of stages required is 5.5. Using the same 
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: Turbine 


No. 1: 


h = 


6 stages 


V. = 
i 


500 ft/sec 


V 

e 


625 ft/sec 


h - 


0.85 


He m 


0.75 


1+f » 


1.00825 


M,/ 


226.5 BTU/lb 


Mi; 


231.8 BTU/lb 


1 T 


0.835 



Since this efficiency is 3,5% larger than assumed in the prelimin- 
ary analysis, a recalculation of exit pressure and specific volume is 
necessary prior to calculating blade heights and diameters at the last 
stage. 




P g = 4 57.6, psia 
V = 8.3 15 ftyib 

A summary of first and last stage data for Turbine No. 2 is listed 
in Table 5-2. 
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V 



TABLE 5-2 





FIRST AND LAST STAGE 


DATA FOR TURBINE NO. 


2 




First Stage 


Last Stage 


Units 


V 2 


1095 


1130 


ft/sec 


V 3 


287 


296 


»t 


°<2 


74.78 


74.78 


degrees 


<*3 


0 


0 


it 


W 2 


408 


421 


ft/sec 


W 3 


825 


850 


m 


@2 


45.29 


45.29 


degrees 


ft 


-67.98 


-67.98 


M 


U 


765 


788 


ft/ sec 


d. 

1 


25.76 


25.76 


in. 


d 

c 


31.88 


32.90 


!t 


d 

o 


38.00 


40.04 


it 


h 


6.12 


7.14 


it 


h/dm 


.192 


.217 




Additional 


Turbine Data: 







N s= 5500 rpm 
7r - 0-854 

* 0.834 

n *s 6 stages 
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CHAPTER VI 



REACTOR 



Symbols used In this chapter: 

Symbol Description 
A total heat transfer area 

CX total coolant flow area 

c fission rate per unit of power 

0 ^ specific heat at constant pressure 

d diameter of coolant channels 

f thermal utilization factor 

g conversion factor , 32.17 



H 

h 

K 2 

k- 

CO 

k 

eff 

L 

L 

o 

M 

m 

N 



or 



height of reactor core 

convection heat transfer coefficient 

size-shape factor 

thermal conductivity 

multiplication constant for reactor 
of infinite size 

multiplication constant for reactor 
of finite size 

thermal diffusion length 

thermal diffusion length for pure material 
mass 

mass rate of flow 

3 

number of nuclei /cm 
number of coolant channels 



Units 

ft 2 

ft 2 

fissions/watt 
BTU/lb °F 
in. 



lb - ft 

m 

lb, 2 

f - sec 

cm or ft 



BTU/hr - ft 

-2 

cm 

BTU/hr - ft 2 



cm 



cm 



grams or lb 
lb/sec 



m 



sec 

°F 

°F/ft 
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Symbol 

P 

o 

P 

Q 

R 

t 

V 

v 

W 

<r 

r 

e 

r 

1 

t 

J 

A 



Description 

power 

pressure 

resonance escape probability 

rate of heat transfer 

radius of reactor core 

temperature 

reactor volume 

specific volume 

coolant velocity 

microscopic cross -section 

absolute viscosity 

fast fission factor 

slowing down area for fast neutrons 

neutrons born per absorption 

flux 

friction factor 
increment of 



Subscript 

a absorption 

av average 

c graphite 

/ 

f fission 

He he 1 ium 

ra mixture of graphite and uranium 



Units 

Watts 

psia 



BTU/hr 
cm or ft 
°F 

3 3 

cm or ft 

ft 3 /ib 

ft/sec 

2 

cm 

lb/hr - ft 

2 

era 



neutrons /cm 



2 



- sec 
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Subscript Description 



u uranium 

w wall 

General 

Although the writer has had limited training and no practical 
experience in the field of reactor technology, an attempt was made in 
this chapter to determine the approximate dimensions of the core of a 
helium-cooled enriched reactor, for the primary purpose of establish- 
ing the pressure drop of the coolant. 

Critical Mass 

In determining the core dimensions and the amount of fuel required 
to produce a critical mass the following assumptions have been made: 

(1) The core consists of a homogeneous mixture of graphite and 

U 235 * 

(2) The core shape is a circular cylinder with a height, H, 
equal to its diameter, 2R. 

(3) The coolant channels occupy 1/3 of the total volume and are 
uniformly distributed throughout the core. 

(4) The power rating of the reactor is based on a cycle thermal 
efficiency of 0.329. 

14 2 

(5) An average flux at full powe* is of the order 10 neutrons/cm 
- sec. 

(6) The effect of reflector and cooling channel liners is not 
included . 

(7) No allowance is made for fuel depletion, poison and tempera- 
ture effects on reactivity. 
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For criticality: 



Keff 



K co Q 



-K 2 r 



+ K a L 2 



(6-1) 



where 



K® — r] 6 p -f 



(6-2) 



K = + $)' 



(6-3) 



= £0-f) 



(6-4) 



f 



i 



I + 



Mrsi 



(6-5) 



Nu fa) u 

To calculate these quantities the following nuclear data is known: 

-24 2 

(Q~a) c =» 0,0045 barns (1 barn =*10 cm ) 

( (Ta) a 650 barns 

v u 

L 2 » 2500 era 2 

o 

6 ~ i.o 

p ~ 1.0 

7 « 2.1 

V « 300 era 2 
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V 



Preliminary calculations indicate that for the power and flux 
requirements : 




14,400 



Therefore : 

f = 0.909 

k„ * 1.91 

2 2 
Lb 227 cm 



By trial and error calculations using equation (6-1): 

K 2 = 0.0013 

H = 159 cm = 5.23 ft 

To allow for the effect of the cooling channels and considering 

that they are distributed uniformly over the entire core, the following 

approximations have been made: 

2 

(1) L increases to 9/4 of its original value 

(2) increases to 9/4 of its original value 

2 

Referring to equation (6-1) the net effect of K is to decrease it 
to 4/9 of its original value or: 

K 2 = 0.000578 

H = 239 cm * 7.85 ft 

The ratio of graphite to uranium by mass is: 



M c 

M, 



1 2 . 

235 



N c 

Na 



zr 736 



Since it can be expected that the density of the graphite will be 



somewhat less than its pure state (1.65 g/cm ), it is assumed to be 
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3 

equal to 1.55 g/cm . Then: 

3 

V . 2/3 ; — * 7.175 X lo' 3 cm 3 = 253 ft 3 

c 4 

M c a 7.175 x 10* 3 x 1.55 * 11.15 x 10 J grams * 24,600 lbs 

M * 11 - 15150 grams * 33.4 lbs 

u 736 



N * 
c 



N 

u 



1.55 x 6.023 x 10 
12 

22 

7,77 x 10 
14,400 



23 



22 

= 7.77 x 10 nuclei/cm 



5.4 x 10^^ nuclei/cm 3 



3 



The power rating of the reactor, based on a 50,000 HP output at a 
thermal efficiency of 0.329, is: 



15 2j OOO HP 
I I 35 X /o 8 Waits 

The average flux is: 



50^00 
O. 32 9 




c P 

N* Vm 



neutrons 

2 

cm - sec 



( 6 - 6 ) 



where C — 3 X |0 



10 



fissions/watt - sec 



01 = 549 



barns 



v. 



rn 



V c — 7/75 X 10 6 cm 3 



or 



p = 1.6 x /o'* 



neutrons 

2 

cm - sec 
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The burn-up rate is: 



U 235 used/da y 



No, of fissions 
day 



No, of absorptions 
fission 



No, of grams 
absorption 



- f) u * 



2 35 

6.023 x / 0* 3 



(6-7) 



= 136 graras/day at full power 

i 

Cooling System 

Although the dimensions of the core and the amount of fuel and 
moderator mixture will increase when allowance is made for the effects 
of fuel depletion, poison and temperature, it is assumed here that this 
will have little effect on the characteristics of the cooling system. 
From the helium properties listed in Table 2-2 for the inlet and exit 
of the reactor and the dimensions of the core: 



li 


i v 


= 126.5' 


ft 3 


a - 


H 


10 

11 


h* 


H 

> 

5° 


m v« v 

a 


- A 9 


ft/sec 



For an average wall temperature of 1300°? the heat transferred to 
the coolant is: 

Q = h A (t w - t J f out - t;„) * 3 600 (6-8) 

a TV/ hr 
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where 



387.8 °F 



t - t SS 

w av 



t -t. = 475.5 °F 

out m 



m « 182.4 lb/sec 



c =1.24 BTU/lb - °F 
P 



or 



,6 



-Li-!®. BTU/ hr - ft -F 
A 



N = 



A 



a ' fa * 



N ^ H 



6090 ' 



d 



r c 



h — (64 d 6 TU/ hr - -f P- 8 P (6-9) 



Another expression for h as a function of d is obtained from 
the turbulent flow heat transfer equation: 



Nu. - 0.02 3 R e aS Ft-° + 



(3-16) 



in which the following gas properties are evaluated for an average 
coolant temperature of 1119°F: 

v = 4.348 ft 3 /lb 

c =1.24 BTU/lb - °F 



/A = 0.10 lb/hr - ft 

k = 0.18 BTU /hr - ft 2 - °F/ft 
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Equation (3-16) reduces to: 



h = ~^| BTU/hr - ft 2 - °F 

d 

From equations (6-9) and (6-10): 



d c 1,08 inches 
h = 178 BTU/hr - ft 2 - °F 

N « 2540 channels 



Pipe Friction 



AP - i 



w 2 H 

2 gVcl 



5 = a/84 Re ' - 



0.022 5 



where R = 36,500 

e 

A P — o. 117 psi 



(6-10) 



(3-13) 



(3-14) 
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TABLE 6-1 



SUMMARY OF THE REACTOR CORE DATA 

H 7.85 ft 

M 15.15 Kg 

u 

M c 11,150 Kg 

v 

V 37? ft 3 

W 49.0 ft/sec 

d 1.08 inches 

N 2,540 channels 

P 0.117 psi (Pipe friction only) 
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CHAPTER VII 



INSTALLATION 



Symbols used in this chapter: 



Symbol 


Description 


Units 


d 


diameter 


in 


g 


conversion factor, 32.17 


lb - 
m 

lb f ' ; 


m 


mass rate of flow 


lb/sec 


P 


pressure 


psia 


t 


thickness 


in 


V 


velocity 


ft/sec 


V 


specific volume 


ft 3 /ib 


A 


increment 


— 


n 


efficiency 




<r 


stress 


psi 


<*> 


defined on page 118 


— 


r 


defined on page 118 


— 


Subscripts 


i 


inside, initial 




o 


outside 




th 


thermal 




w 


wall 
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General 



In the previous chapters each of the major components of the cycle 
has been analyzed. The following data were obtained: 

(1) Basic internal dimensions of the turbo-machines such as 
diameters and blade heights. 

(2) Turbo-machine characteristics such as stage and overall 
efficiencies, number of stages, and velocity diagrams at the 
mean radius for the first and last stages, 

(3) Heat exchanger dimensions including dimensions and arrange- 

H 

ment of tubes, shell inside diameter, and total internal 
volume . 

(4) Heat exchanger pressure losses due to pipe friction. 

In order to combine these components in an installation, addition- 
al information is necessary. Such items as piping, casings, inlet and 
exit connections, bearings, couplings, auxiliary equipment, control 
apparatus, and expansion joints, must be designed. It is beyond the 
scope of this investigation to cover all of these itejns, so only those 
which affect the size of the plant and pressure drops will be consider- 
ed in this chapter. 

Before considering the individual components there follows a 
review of some of the relations used in this chapter. 

* * 

Effect of Changes of Specific Volume between Inlet and Exit of 
Heat Exchangers 

In the calculations for pipe friction using equation (3-13) a mean 
specific volume was used. Given the channel below: 
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////// ///////////.//////// 

2 ^ / 

7777777777777 777777777777 



Figure 7-1 Channel of Constant Flow Area 



A 

the pressure change between (1) and (2) due to the change in specific 
volume can be approximated by the equation: 



R - P, 



(vf M~ V 2 

\ V /av. Q x 1 4 4 




(7-1) 



Pressure Drop for Flow Acceleration 



Vi 



\\\ \ 



\\ \\\\\\ 

\ 



Vi. 








Figure 7-2 Channels of Decreasing Flow Area 





\ 2 < 3 V x 144 



psi 



(7-2) 
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where (p is the efficiency of conversion from static 
pressure to velocity head. 



Pressure Increase for Flow Deceleration Due to a Sudden Enlargement 

.WWW .W 










t 



\ A \ \ \ \ \ 



Figure 7-3 Channel with an Abrupt Increase in Flow Area 



From the momentum theorem, 



PrP, 



v, - v/ 

9 v * I 44 



pSL 



Pressure Increase in a Diffuser 



WW\\ 






\\ \ \ 



M - 



\- WAW 



VX'VTW 



Figure 7-4 Diffuser 



(7-3) 






y V/) 



P5L 



(7-4) 



2gv V H4 

/ 

where is the efficiency of conversion from velocity 
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head to static pressure. 



Pipe Dimensions 

The internal pipe diameters and flow velocities are related as 
follows : 



di = (p-^r) 2 * 12 lnches 



(7-5) 



The wall thickness, t , depends on the allowable circumferential 
stress, : 



t 



w 



(a P) dj 

Z<r 



inches 



(7-6) 



y jpAgg 

In analyzing the turbines and compressors inlet and exit pipe 

* 

velocities were specified in order that overall efficiencies could be 
determined. The internal pipe diameters corresponding to these velo- 
cities can be calculated . from equation (7-5). It is desirable to use 
standard pipe sizes although these may differ slightly from the calcu- 
lated diameters. Using the pipe property chart and allowable stress 
data from Littleton^^ , Table 7-1 was compiled with the use of equa- 
tions (7-5) and (7-6). 

The actual outside diameters of some of the installed piping will 
be greater than those listed by the amount of insulation required. For 
locations four through eight (Figure 2-1) this will be in the order of 
six to eight inches. For the other locations none will be required 
unless some is desired to limit the space temperature. 
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I'll'lNi: DATA 




I In fit »m v< 1 1 >» l I 1 on n a itimiotl tin I ho i y« I n /iml il t i I n i n l t nh( I y t i out I limit k o i i » wi poml t In 

1 liip | ' I pt* <1 i mm’ t i i m I i « I ml . 



Regenerator and Precooler 



It is convenient to discuss these two units together since they are 
in series for cooling the hot gas between the power turbine discharge 
and the inlet to Compressor No. 1. Figure 7-5 illustrates the two units 
showing approximate dimensions. For an allowable stress of 10,000 psi 
for the regenerator shell, its thickness is two inches. An allowable 
stress of 12,000 psi for the precooler shell requires a thickness of 
one inch. For purposes of calculating pressure changes, schematics of 
the flow passages are illustrated in Figures 7-6 and 7-7. 

Referring to Figure 7-6, the flow has been broken down into the 
► 

following parts; 

(a) Deceleration by means of a diffuser of 2 foot length between 
(1) and (2), with corresponding pressure recovery. 

(b) Deceleration to zero velocity at (3) with no pressure recovery. 

(c) Acceleration to V^. 

(d) Pipe friction effects between (4) and (5). 

(e) Deceleration from to zero velocity at (6) with no pressure 
recovery. 

(f) Acceleration to V-,. 

(g) ~ (k) For remainder of the channel between (7) and (12) the 
effects are similar to those outlined in items (b) - (f) above. 

(1) Effects of changes in specific volume between inlet and exit. 

A summary of the velocities and pressure changes for the locations 
specified in Figure 7-6 is listed in Table 7-2. Where applicable, 

V" = 0.85 and (p = 0.80. 
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Figure 7-5 Regenerator and Precooler 




5 



4 



3 



6 




Figure 7-6 Schematic Representation of Flow through the 

Regenerator and Precooler 





Figure 7-7 Schematic Representation of Flow through the 

Regenerator (H.P. Side) 
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TABLE 7-2 



SUMMARY OF 


PRESSURE 


LOSSES IN THE 


REGENERATOR AND 


P RECOOLER 


Location (x) „ 

(Figure 7-6) v(ft J /lb) 


V( ft/sec) 


P - P ,(psi) 
x x-1 r 


Equation 


1 


8.413 


625 







2 


ft 


458 


1.97 


(7-4) 


3 


)l 


0 


0 


(7-3) 


4 


tf 


53.6 


-0.05 


(7-2) 


5 


4.607 


28.7 


- 0.76 


(3-13) 


6 


tl 


0 


0 


(7-3) 


7 


It 


200 


-1.17 


(7-2) 


8 


It 


0 


0 


(7-3) 


9 


tt 


86 


-0.22 


(7-2) 


10 


3.428 


64 


-0.72 


(3-13) 


11 


tt 


0 


0 


(7-3) 


12 


It 


400 


-6.30 


(7-2) 


Comp ressibility 


Effect 




0.12 


(7-1) 






A "total 


= -7.13 








A p 

P. 

1 


= 0.0158 
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TABLE 7-3 



SUMMARY OF PRESSURE LOSSES IN THE REGENERATOR (H. P. SIDE) 

Location (x) P - P fnsO 

(Figure 7-7) v(ft /lb) V(ft/sec) x x-P p w Equation 



1 


1.883 


450 


— 


— 


2 


11 


323 


7.22 


(7-4) 


3 


II 


0 


0 


(7-3) 


4 


II 


66 


-0.31 


(7-2) 


5 


3.677 


132 


-12.03 


(3-13) 


6 


If 


0 


0 


(7-3) 


7 


II 


600 


-13.20 


(7-2) 


Compressibility Effect 




-1.09 


(7-1) 






A P 


= -19.43 








AP 

P. 

t 


« 0.0194 




Considering 


the conditions inside the 


tubes of the 


regenerator, 


the schematic of 


Figure 7-7 


indicates that 


items (a) - ( 


f) on page 121 



apply. The results are listed in Table 7-3. Where applicable, 
Y' » 0.85 and (f) * 0.80. 
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Figure 7-9 




Figure 7-8 Intercooler 




Schematic Representation of Flow through the Intercooler 
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Intercooler 



The intercooler is similar to the precooler and the flow is broken 



down in the. same 


manner . 


The dimensions 


indicated in Figure 


7-8 are 


approximate and 


based on 


calculations in 


Chapter III. For an 


allowable 


stress of 14,000 


psi, the 


shell thickness 


is one inch. 








TABLE 7-4 






SUMMARY OF PRESSURE LOSSES IN 


THE INTERCOOLER 




Location (x) 
(Figure 7-9) 


v(ft 3 /lb) V(ft/sec) 


P K - F X-l' P8l > 


Equation 


1 


2 . 768 


400 


— 




2 


II 


294 


2.43 


(7-4) 


3 


II 


0 


0 


(7-3) 


4 


II 


82 


-0.32 


(7-2) 


5 


2.309 


68 


-1.05 


(3-13) 


6 


II 


0 


0 


(7-3) 


7 


II 


300 


-5.26 


(7-2) 


Compressibility 


Effect 




0.09 


(7-1) 






A P 


= -4.11 








A P 

P. 


= 0.0062 
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Reactor 



As illustrated in Figures 7-10 and 7-11, the flow through the 
reactor is similar to that through the heat exchangers. The overall 
dimensions are only estimates, since no calculations were made concern- 
ing the thickness of shielding and reflector required. 

TABLE 7-5 

SUMMARY OF PRESSURE DROPS IN THE REACTOR 



Location (x) 


v(ft 3 /lb) 


V(f t/sec) 


P - P .(psi) 
X x-1 


Equation 


i 


3.677 


600 










2 


fl 


434 


4.29 


(7-4) 


3 


1! 


13.8 


0.34 


(7-3) 


4 


tl 


41.5 


-0.06 


(7-2) 


5 


5.019 


56.5 


-0.12 


(3-13) 


6 


II 


18.9 


0.03 


(7-3) 


7 


It 


600 


-9.69 


(7-2) 


Compressibility 


Effect 




-0.04 


(7-1) 








= -5.25 








Ap 

p 


= 0.0054 








L , 
1 
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Figure 7-10 Reactor Figure 7-11 Schematic Representation of 

Flow through the Reactor 




Tables 7-6 and 7-7 compare the calculated machine efficiencies 



and pressure losses with those assumed in Chapter II. Another cycle 
analysis using these calculated values results in the following new 
performance indices : 



m 

th 

helium rate 
work ratio 



162.5 lb/sec (For 50,000 HP output) 
0.361 

11.7 lb/HP - hr 
0.424 



TABLE 7-6 



COMPARISON OF ASSUMED AND CALCULATED VALUES OF AP/P. 



Unit 



AP/P i (assumed) 



A P/P^ (calculated) 



Regenerator 



High pressure 
side 


0.015 


0.0194 


Low pressure 
side 


0.010 


0.0158 


Frecooler 


0.0075 




Reactor 


0.010 


0.0054 


Intercooler 


0.005 


0.0062 
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TABLE 7-7 



COMPARISON OF ASSUMED AND CALCULATED VALUES OF MACHINE EFFICIENCIES 



Unit 


ij (assumed) 


Ij (calculated) 


Compressor No. 1 


0.84 


0.867 


Compressor No. 2 


0.80 


0.840 


Turbine No. 1 


0.85 


0.858 


Turbine No. 2 


0.80 


0.835 



n 




s 



M 

O 

4J 

a 

c4 



Figure 7-12 Power Plant Arrangement 





CHAPTER VIII 



STRESSES 



Symbols 


used in this chapter: 


Symbol 


Description Units 


A 


. 2 

area in 


b 


width of turbine disk in 


c 


chord in 


c 

P 


specific heat at constant pressure BTU/lb - °F 


F 


force lb^ 


S 


lb 

conversion factor, 32.17 m ft 

lb, 2 

f sec 


h 


blade height in 

or 

enthalpy BTU/lb 


I 


4 

moment of inertia in 


J 


conversion factor, 778 ft - lb/BTU 


K 


coefficient defined in equation (5-10b) 


M 


bending moment in - lb 


m 


mass lb 

m 


m 


mass rate of flow Ib/sec 


N 


rotational speed rpm 


P 


pressure psia 


q 


centrifugal loading lb/in 


r 


radius in 


s 


spacing in 


T 


absolute temperature °R 
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Symbol 



t 

U 

V 

v 

W 

Z 

oc 

£ 

Y y' 

5 , 6 ' 

<p 

<T 

co 

A 

(° 

X 

7 



Description 

temperature 

peripheral speed 

absolute velocity 

specific volume 

relative velocity 

no. of blades 

absolute flow angle 

relative flow angle 

angles defined by Figure 8-3 

angles defined by Figure 8-2 

* — 
acute angle between y-axis and F 

acute angle between F^ and m 

stress 

rotational speed 

increment 

density 

degree of reaction 
efficiency 



Subscrip t 

c 

CG 

g 

i 

m 

min 



Definition 

centrifugal or center 

center of mass 

gas force 

inner 

axial 

minimum 



Units 

°F 

ft /sec 
ft /sec 
ft 3 /ib 
ft/sec 

degrees 

degrees 

degrees 

degrees 

degrees 

degrees 

psi 

rad/sec 

lb/ft 3 
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h 






Subscripts Definition 



max 



maximum 



o 



outer 



R 



rotor 



s 



stator, isentropic, stage 



x, y 



coordinates 



overscore vector 

General 

Of the quantities that were assumed in the analysis of the turbo- 
machines the most critical were the allowable peripheral speeds. These 
had an effect on the machine efficiencies, radial dimensions and the 
number of stages required. The choice of the allowable speeds was made 
for both the compressors and turbines on the premise that the blade and 
rotor stresses would not exceed safe limits. It is the objective in 
this chapter to investigate the parts of the machines where the stresses 
will be the greatest; namely, the first and last stages of Compressor 
No, 1 and the first stage of Turbine No. 1. 

Compressor Blades 





Figure 8-1 Forces Acting on a Compressor Rotor Blade 
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For free-vortex flow and equal work input at all radii: 



rV « constant 

u 

V = constant 

m 

From the first and last stage data given in Tables 4-2 and 4-5, 
the equations above, and Figure 4-3, the velocity diagrams for the 
inner, center and outer radii have been determined as indicated in 
Tables 8-1 and 8-2. 
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TABLE 8-1 



VELOCITY 


DIAGRAM DATA 


FOR THE FIRST 


STAGE OF 


COMPRESSOR NO. 




Inner 


Center 


Outer 


Units 


r 


19.435 


20.25 


21.065 


in . 


v 1 *v..*=v 

13m 


436 


436 


436 


f t/se 


V 2 


616 


604 


593 


II 


V 

U 1 


0 


0 


0 


II 


V 

U 2 


435 


418 


402 


It 


W 1 


1263 


1310 


1358 


(1 


W 2 


868 


928 


987 


II 


w 

U 1 


1185 


1235 


1285 


II 


w 

U 2 


750 


819 * 


883 


II 


u 


1185 


1235 


1285 


IS 


Z 3 ! 


69.8 


70.6 


71.2 


degre 


@ 2 


59.8 


62.0 


63.7 


it 


04 1 


0 


0 


0 


n 


Od 

2 


45.0 


43.8 


42.7 


it 


c 


0.82 


0.82 


0.82 


in. 


S R 


0.79 


0.82 


0.85 


n 


s 


0.37 


0.39 


0.41 


ii 



s 
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TABLE 8-2 



VELOCITY DIAGRAM DATA FOR THE LAST STAGE OF COMPRESSOR NO. I 





Inner 


Center 


Outer 


Units 


r 


19.435 


20.10 


20.765 


in . 


V,=V,=V 
I 3 m 


432 


432 


432 


ft/sec 


v. 


609 


599 


591 


»i 



\ 0 0 0 



V 

u 2 


429 


415 


402 


ft 


W 1 


1262 


1300 


. 1333 


tl 


W 2 


869 


920 


964 


It 


W 

U 1 


1185 


1225 


1267 


tl 


w 

U 2 


756 


810 


865 


If 


U 


1185 


1225 


1267 


ft 

V 


Pi 


70 


70.6 


71.1 


degrees 


CNI 


60.2 


62.0 


63.4 


it 


**1 


0 


0 


0 


>* ii 


U 2 


44.8 


43.8 


43.0 


h 


c 


0.82 


0.82 


0.82 


in. 


S R 


0.79 


0.82 


0.85 


ti 


s 

s 


0.38 


0.39 


0.40 


it 



Blade Profiles 

The profiles selected for the rotor and stator blades are NACA 
Airfoil 6309 and modified Joukowsky Airfoil 603, respectively. These 
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are illustrated in Figures 8-2 and 8-3 showing important dimensions. 
Table 8-3 lists important section properties. Angles T and 8 
depend on space/chord ratio, the velocity diagrams, and the type of 
profile. Actual values of 8 were obtained from de Haller whereas 
values for ^ were obtained from Vavra^^. These are tabulated in 
Table 8-4. 

TABLE 8-3 

AIRFOIL SECTION PROPERTIES FOR A CHORD OF 0.82 INCHES 









NACA 6309 




Joukousky 603 






Area 




0.042 




0.0532 


. 2 
m 




I . « 

mm 


I 

>: 


1.46 x 10 5 




5.4 x 10" 5 


. 4 
in 




I 

max 


I 

y 


1.565 x lO -3 




1.89 x 10" 3 


in 4 








TABLE 8-4 












VALUES OF f 


AND 6 (FIGURES 8- 


•2 AND 8-3) 








First Stas 


;e. 




Last Stage 






Inner 


Center 


Outer 


Inner 


Center 


Outer 


6 


27.2 


25.2 


23.5 


27.3 


25.3 


23.8 


6' 


27.5 


25.5 


23.8 


27.6 


25.6 


24.1 


r 


26.0 


25.0 


24.8 


— 


— 


— 


'T 


/ 

29.1 


28.1 


27.9 
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NACA 6309 



1 - 




Location 


x (in) 


y (in) 


A 


0.0928 


0.0492 


B 


0.347 


-0.041 


C 


0.473 


-0.034 



Figure 8-2 First Stage Rotor Blade for Compressor No. 1 



140 




Figure 8-3 First Stage Stator Blade for Compressor No. 1 
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Centrifugal Blade Forces 



The maximum blade stress will occur at the blade root. A portion 
of this stress for the rotor blade will be due to the centrifugal force 
effects and can be calculated from the following equations: 

F c - to 2 l hi- (8-1) 

07 = 4 - ( 8 - 2 ) 

A 

3 

For constant chord blades of low alloy steel ( * 0.282 lb/in ) 
the centrifugal forces and stresses are as listed in Table 8-5. 



TABLE 8-5 

CENTRIFUGAL FORCES AND CORRESPONDING STRESSES FOR 
THE ROTOR BLADES OF COMPRESSOR NO. 1 





First Stage 


Last Stage 


Units 




733 


733 


rad/sec 


r CG 


20.25 


20.10 


in. 


m 


0.0193 


0.0158 


lb 

m 


F 

c 


544 


442 


lb f 


<Tc 


12,950 


10,500 


psi 


Gas 


Bending Force 







The force acting on the blades due to the pressure differential 
and the change in momentum of the flow can be calculated by means of 
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I 




the momentum theorem. 




Figure 8-4 Momentum Theorem Control Surfaces for the Blades 

of Compressor No. 1 



Referring to Figure 8-4 and using the dashed line as a control 
surface, the momentum equations become: 

For Rotor: 

- If* W ' + 'T* = A . _ - F 3 (a-3) 

where A, = A 2 = Sh 

For Stator: 



^z 

3 









3 

where /\ „ 



V, = (HA* - P.A^in-F, 

= A, = sh 



(8-3a) 
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blades, 2T > and the pressure rise across the blades in addition to 
quantities previously determined. 



The pressure rise for the first stage is easily calculated from 
equation (l-6a) since and are known from Table 2-2 and A h g from 
Table 4-4. Although the degree of reaction, pc, is defined in terras of 
enthalpies, it can be expressed in terms of pressure rise with little 
error . 



For the last stage, the quantities A h and T^ must be calculated. 
This can be done by equations (4-24) and (l-5a) since h^, U, T* and Yj s 
are known. The pressure rise can then be calculated from equations 
(l-6a) and (4-15). 




(8-4) 



;< 



hz ~ hj 

h 3 -h 




(4-15) 



Table 8-6 indicates values of Z , 



AP and F . 

g 
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TABLE 8-6 



GAS FORCES ACTING ON THE 
First Stage 





Rotor 


Stator 


z 


155 


326 


A P 


24.2 


4.9 


F 

g 


36.0 


7.92 


e 


26.2 


66.7 




The gas bending 


stresses have 


on 


the rotor blades 


and one on the 



Referring to Figures 8-2 and 

<r = _ Mil 

9 T T 



BLADES OF COMPRESSOR NO. 1 



Last Stage 


Units 


Rotor 


Stator 




144 


324 


Blades 


29.8 


6.1 


psi 


35.9 


lbs 


25.2 


degrees 



been calculated for three locations 
first stage stator blade. 
t-3, the stresses can be expressed as: 

(8-5) 



where = — y S in (J> 

= - § Pcj Cos $ 

Q — $ {or rotor 

~ 90 ) t y' for stator 
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TABLE 8-7 



BLADE ROOT STRESSES DUE TO GAS FORCES 





Rotor 




Stator 






First Stage 


Last Stage 


First Stage 


Units 


$ 


-1.3 


-2.4 


-*4.6 


degrees 


M 

X 


-29.3 


-23.8 


-6.42 


in - lb 


M 

y 


i 0.62 


1.0 


+0.52 


in - lb 


( Va 


-99,000 


-80,300 


+12,850 


psi 


<Vb 


+82,350 


+66,600 


— 


psi 


(ti Vc 


+68,500 


+55,800 


— 


psi 



Applying the centrifugal force stresses to the gas force stresses, 
a maximum stress of 95,300 psi occurs at point B of the first stage 
rotor blade. This stress is higher than normal long-life design criter- 
ia and would preferably be reduced to 50,000 psi or less. Lower 
stresses can be obtained by modifying the blades. Instead of using a 
9% thick blade at all sections, a tapered blade 12% thick at the root 
would reduce the stress to 53,700 psi. If, in addition, the chord is 
increased by 10%, the stress is further reduced to about 40,400 psi. 

There are many low-alloy, high-strength steels that would be suit- 
able for the compressor blades such as AISI Nos. 4140, 4340, 6150, 

8640 and 8740. 

As one might expect, the stresses in the stator are small compared 
to the rotor stresses since most of the pressure rise is in the rotor 
blades. Providing there is no serious vibration problem, no major 
design changes would have to be made to permit a mean peripheral speed 
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of 1235 ft/sec. 



Turbine Blades 




Figure 8-5 Forces Acting on a Turbine Rotor Blade 



For free-vortex flow and equal work output at all radii: 



rV * constant 
u 



V «= constant 
m 



From the first stage data given in Table 5-1, the equations above, 



and Figure 5-3, the velocity diagrams for the inner, center and outer 



radii have been determined as indicated in Table 8-8. 
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TABLE 8-8 



VELOCITY DIAGRAM DATA FOR FIRST STAGE OF TURBINE NO. I 





Inner 


Center 


Outer 


Units 


r 


10.945 


13.575 


16.205 


in. 


V 2 


1370 


1121 


9.52 


ft/sec 


W V m 


294 


294 


294 


M 


V 

U 2 


1339 


1080 


904 


ii 


V 

U 3 


0 


0 


0 


n 


W 2 


763 


418 


296 


it 


W 3 


696 


847 


982 


n 


W 

u 2 


706 


295 


-33 


n 


W 

U 3 


-633 


-785 


-937 


n 


U 


633 


785 


937 


it 


/3 2 


67.4 


45.29 


-6.4 


degrees 




-65.0 


-67.98 


-72.6 


h 


°<2 


77.6 


74.78 


72.0 


it 




0 


0 


0 


ti 


S R 


0.81 


1.0 


1.19 


in. 


9 

S 


0.97 


1.2 


1.43 


it 



Also: N = 7000 rpm 

h = 5.26 in. 

Z R - 85 

Z - 71 

s 
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Rotor Blades 



The root, center and tip profiles of the first-stage rotor blade 
are illustrated in Figure 8-6. This is a tapered blade in which the 
centroids of each section lie on a radial line. The maximum blade 
stresses occur at the blade root, hence only these section properties 
have been determined. As in the compressor, the stress is due to a 
combination of gas bending force and centrifugal force effects* 
Centrifugal Force 

The volume and center of mass of the tapered blade having root, 
center and tip sections as shown in Figure 8-6 can be determined from 
the following plot: 




Blade Volume 



0.21 + 0.48 , _ r 3 

- x 5.26 *s 1.82 in 



CG 



10.94 + 2.30 » 13.24 in 



For alloy S-816: ( fO = 0.310 lb/in ) 



m = 0.564 lb per blade 

F c « 10,500 lb per blade 

G~ * 21,900 psi 
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Tip 



Center 



Root 




Figure 8-6 First Stage Rotor Blade Profiles for Turbine No. 1 



150 



Gas Force 



Before the gas force can be calculated by mean9 of the momentum 
theorem, the pressure differential between the inlet and exit to the 
rotor blade must be calculated. 

Referring to Figures 5-4 and 5-5, the following relations are 
apparent: 






v; 



V/ 



(ff 2gJ 



h 2 + 



2\ w; 



J 



b, ~ 



2g J 



(Ah;)" — (Ah s y 



Vp' 



<T 



Also ; 







/pOft 1 




= T 


- (1 J 



(8-6) 



(5-12) 



(8-7) 



(8-8) 



(l-7a) 



V = 



- h 



6oo 



/. 2f 



+ 600 



(1-5) 
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(5 -10a) 



/ , N * Ah s 

/ S 



(Ah s ) — (Ah^ + K 



V' 






(5-10b) 



Combining equations (8-8), (5-10a) and (5-10b): 



(Ah y^-k 

I s ' ^ ^ 2g<T ^ 2g-J 



(8-9) 



Also : 



W = 7 c, 






pf 



(l-7a) 



From these equations, Tables 2-2, 8-8 and Chapter V, the pressure 
rises are: 

?1 - P 2 = 27,7 psi 

P 2 - ? 3 a 11.5 PSI 

Using the dashed line of Figure 8-7 as a control surface, the 
momentum theorem becomes : 






k + -4 1 w 3 = (ev6A>-e 



where A n » 



= sh 



or 



F i 

F 9 



— 60.5 m +- 72.0 U 

- 94-. 3 lb • 9 = 49.9 
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Figure 8-7 Momentum Theorem Control Figure 8-8 Turbine Nozzle Profile 

Surface for a Turbine 
Rotor Blade 



The fibers most heavily stressed in the root cross section are 
indicated by A and B in Figure 8-6. Axes x and y are principal axes. 



_ AVj _ M,* 

T ' 

where ^ E (°S 4 1 

-k- Rj Sin 4 



(p = 72.5° - © => 22.6° 

— 70 30 p5L 

(fy s = 66/0 P 5C< 

A total stress of 29,000 psi exists at A. Jennings states that 
the blade root temperature will normally be 175-300 a F lower than the 
gas stagnation temperature or a temperature of about 1100°F here. The 
stress -to-rupture in 100,000 hours for S-816 alloy is 47,000 psi at 
1100°F. Therefore, the rotor blade stresses are not excessive for the 
blading of Figure 8-6 providing no serious vibration problem exists. 
Stator Blades 

There is no attempt here to do a detailed stress analysis on the 

nozzles. Figure 3-8 illustrates a typical profile for the root section. 

The gas force has been calculated to be 200 lbs at 0 a 25.6°. The 

minimum section modulus, *nin , of the profile is approximately 17.6 x 

y 

max 

•3 3 

10 in and can be increased by increasing the % thickness and/or 
chord length. 
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30,000 psi 




I* 



Since the stator blade temperatures will be higher than the rotor 
blades, a maximum stress of less than 25,000 psi is desirable and this 
can easily be achieved by slight modification of the blade. 



Turbine Disks 





(«) 



(b) 



Figure 5-9 Turbine Disk 
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The type of disk being considered is illustrated in Figure 8-9. 

All dimensions except the disk width, b, are known from blade data and 
practical considerations. The value of b depends on stress conditions 
due to the centrifugal forces. The stresses in the rim (r A <r<r.) are 
substantially smaller than the blade root stresses so they need not be 
calculated here. Since we are interested in the stresses and dimensions 
of the remainder of the disk, the blade and disk assembly can be replac- 
ed by the hub portion only (0<r<r^) with an external load due to 
centrifugal forces of the mass of the blades and the rim (Figure 8-9(b)). 

Because of the irregular shape of the rim it is divided into sec- 
tions (Figure 8-10) for the purpose of determining centrifugal forces. 




r = 8 . 695 " 

A 



Figure 8-10 Approximate Cross Section of Turbine Disk Rim 



TABLE 8-9 



CENTRIFUGAL LOADING ON THE TURBINE DISKS 



Item 


Material 


P lib /in 3 ) 


r CG (in) 


ra(lb/in at A) 


q(lb/in 


i 


Fir Tree: S-816 
Rim: 16-25-6 


( 

0.300 


10.145 


0.957 


13,600 


2 


16-25-6 


0.291 


9.170 


0.150 


1,920 


3 


16-25-6 


0.291 


8.845 


0.266 


3,270 


4 


16-25-6 


0.291 


8.995 


0.126 


1,580 




Loading at A due to 


blades : 




rt * 

•rim 


20,370 




q(blades ) 


F Z 
= c 


« 16,350 


lb/in at A 





Total loading: 

q(total) = 16,350 + 20,370 » 36,720 lb/in at A 



For a constant stress disk: 



lr A z -r^) 

b = 



(8-10) 



Since the disk temperature will be approximately 1000°F at A and 
500°F at the center, an allowable stress of 25,000 psi is assumed, 
based on 100,000 hour stress-to-rupture data of 16-25-6 alloy. 




3 6 720 
25000 



— 1.47 i n . 



at A) 
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Values of b for various radii are tabulated in Table 8-10. 



TABLE 8-10 



WIDTHS OF 


CONSTANT STRESS DISK 


r ( in -) 


b (in.) 


0 


3.15 


1.0 


3.12 


2.0 


3.02 


3.0 


2.87 


4.0 


2.68 


5.0 


2.45 


6.0 


2.19 


7.0 


1.96 


8.0 


1.65 


8.695 


1.47 



Figure 8-11 illustrates the first and second stage rotor and 
nozzle assemblies. 

It appears from the foregoing stress calculations that the peri- 
pheral speed that was assumed is not excessive. There is, of course, 
the possibility of much higher blade stresses if the unsteady gas force 
frequency is close to the blade natural frequencies. This problem will 
not be investigated here since normally these effects can be minimized 
by modification of the blading rather than major redesign of the 
turbine . 
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We/c/ 




Figure 8-11 First and Second Stages of Turbine No. 1 
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CHAPTER IX 



* 



SUMMARY AND CONCLUSIONS 



In this report, a limited design study was performed for each of 
the major components of a 50,000 HP closed-cycle gas turbine power 
plant employing helium a3 the working fluid and a nuclear reactor as 
the heat source. The descriptive data tabulated in Table 9-1 should be 
considered representative rather than optimum for a cycle having a maxi- 
mum pressure of 1000 psia and a maximum temperature of 1350°F. It is 
believed that comparable sizes and performance would be obtained if 
different choices are made for such variables as compressor and turbine 
rpm, type of blading, ratio of blade spacing to chord, type of heat 
exchanger surfaces, and reactor cooling channel arrangement. 

The thermal efficiency of the cycle is in the order of 33-367* 
depending on the accuracy of the machine efficiencies and pressure 
drops. For temperatures of 1450 - 1500°F it is estimated that the 
cycle efficiency would increase to about 407>. The size of the compon- 
ents would be reduced substantially from those herein if higher pres- 
sures were employed. 

Although this design study was based on a ship propulsion plant, 
no mention has been made as to its operational requirements. One of 
these requirements is astern power. It may be possible to use an astern 
gas turbine in an arrangement similar to that of a conventional steam 
propulsion plant. 

Another operational problem of critical importance concerns the 
control of the cycle. The primary control consists of raising and 
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TABLE 9-1 



SUMMARY OF DESCRIPTIVE DATA FOR THE MAJOR COMPONENTS 



Unit 



Stages RPM 



Max. Blade Tip 
Diameter (in.) 



Max. Blade 

Height (in.) Efficiency 



Compressor #1 


7 


7000 


42.1 




1.6 


0.867 


Compressor #2 


10 


7000 


40.9 




1.5 


0.840 


Turbine #1 


9 


7000 


34.3 




6.2 


0.858 


Turbine #2 


6 


5500 


40.0 




7.1 


0.834 




Overall 

Length 

(ft) 


Outer 

Diameter 

(ft) 


Tube 

Length 

(ft) 


Tube 
I. D. 
(in. ) 


No . of 
Tubes 


Ail 

b. 

l 


Regenerator 


45 


7.5 


35.7 


0.30 


10,330 


0.0194 


Precoolcr 


18 


4.6 


10.2 


0.40 


2,690 


0.0158* 


Intercooler 


18 


3.8 


9.9 


0.40 


1,700 


0.0062 


Reactor Core 


7.85 


7.85 


7.85 


1.08 


2,540 


0.0054 



* Includes the low pressure side of the regenerator. 



lowering the upper pressure level by means of an auxiliary system con- 
taining high pressure and low pressure helium accumulators. The speed 
with which this control functions will have an effect on ship maneuver- 
ability. Due to the large inventory of helium in the cycle this method 
of control may be sluggish compared to that of conventional steam 
plants. One important characteristic of this type of cycle is its good 
efficiency over a wide range of power. 

Since helium is unavailable at sea except for a limited reserve 
supply, there must be practically zero leakage. This is very difficult 
to achieve when rotating machinery is involved. The main leakage prob- 
lem occurs at Turbine No. 2 where the shaft leaves the casing. The 
remainder of the cycle can be enclosed in a leakproof casing. Bearing 
lubrication then becomes a problem. If oil is used as the lubricant, a 
rather complex system of separating helium from the oil is required in 
order that the helium can be returned for further use. It is possible 
to design the bearings so that high pressure helium is the lubricant 
while the machines are operating. With this arrangement, an auxiliary 
oil system would be required for starting. 

Long half-life fission products must be prevented from being car- 
ried to all parts of the cycle by proper design of the cooling channels 
and filtering system. Then, only the reactor would have to be shield- 
ed. 

These are only some of the problems that would have to be thor- 
oughly investigated before a meaningful appraisal of the closed-cycle 
helium gas turbine plant can be made. However, on the basis of the 
limited investigation of this report, such a cycle has definite advan- 
tages and is worthy of further study. 
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